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1.  INTRODUCTION 


Air  Force  Contract  F19628-90-K-0026  was  established  June  29, 1990  at  Louisiana 
State  University  to  support  work  leading  to  a  model  of  the  heavy  ion  environment 
encountered  by  the  CRRES  (Combined  Release  and  Radiation  Effects  Satellite) 
spacecraft.  Such  a  model  is  needed  to  support  the  Microelectronics  Package  data 
analysis  and  to  understand,  quantitatively,  the  effects  of  the  near-Earth  space 
environment  on  electronic  devices.  The  overall  effort  involves  work  by  three 
laboratories;  The  University  of  Chicago  (J.  A.  Simpson  and  M.  Garcia-Munoz,  lead 
investigators).  The  Lockheed  Palo  Alto  Research  Laboratory  (D.  L.  Chenette,  lead 
investigator)  and  Louisiana  State  University  (J.  P.  Wefel  and  T.  G.  Guzik,  lead 
investigators).  The  principal  source  of  data  on  high  energy,  heavy  ions  on-board 
CRRES  is  The  University  of  Chicago  Experiment,  ONR-604,  (J.  A.  Simpson,  PI;  M. 
Garcia-Munoz  and  J.  P.  Wefel,  Co-I’s)  which  was  designed  specifically  to  measure 
the  charge,  mass  and  energy  of  particles  from  Hydrogen  through  Nickel  that  come  to 
rest  in  the  instrument.  Other  sensors  on  CRRES  provide  detailed  measurements  of 
protons  and  electrons. 

The  overall  goal  of  this  effort  is  to  provide  a  numerical  model  of  the 
high  energy  (>10  MeV/nucleon)  interplanetary  heavy  ion  environment, 
which  is  operationally  interchangeable  with  the  existing,  widely-used 
computer  model  CREME  (see  Adams  et  al.,  1981;  1982).  This  new  model  will 
be  updated  periodically  using  data  from  the  charged  particle  sensors  on  the 
CRRES  satellite.  The  effort  involves  (a)  modifications  to  existing  computer 
codes  to  include  time  dependent  parameters  that  can  be  adjusted  to  reproduce 
the  measured  data,  (b)  timely  processing,  reduction  and  analysis  of  the  data, 
including  uncertainty  estimates,  from  selected  CRRES  sensors,  (c)  comparison 
of  the  experimental  results  to  the  model  predictions  and  adjustment  of  the 
model  parameters  and  (d)  preparation  and  dissemination  of  periodic  reports 
giving  the  analysis  results  and  the  refined  model  parameters  to  the  user 
community. 

This  approach  to  the  assessment  of  the  heavy  ion  environment  provides  both 
consistency  in  technique  for  the  user  groups,  especially  during  the  analyses  of  the 
CRRES  Microelectronics  data,  and  current,  accurate  information  on  the  heavy  ion 
environment  deduced  from  instruments  on-board  the  same  spacecraft,  which  will 
be  consistent  with  the  accumulated  body  of  experimental  data  on  interplanetary 
heavy  ions. 

The  project,  "The  Interplanetary  Heavy  Ion  Environment  for  both  Galactic 
Cosmic  Rays  and  Solar  Energetic  Particles  during  the  CRRES  Mission,"  is  a  multi¬ 
year  effort.  This  Technical  Progress  Report  covers  the  first  15  months  of  effort, 

~1  July  90  *  30  September  91.  The  first  model  release  is  now  scheduled  for  Fall,  1991 
in  the  mid-November  to  early  December  time-frame. 
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2  THE  "HEAVY  ION  MODEL"  PROJECT 


a.  Rationale 

* 

The  heavy  ion  radiation  environment  encountered  in  near-Earth  space 
consists  of  three  major  components:  (a)  galactic  cosmic  rays  (GCR),  (b)  solar 
energetic  particles  (SEP),  and  (c)  trapped  magnetospheric  particles.  The  first  two  of  » 

these  arrive  at  the  earth  from  outside  our  immediate  geospace  environment  and  are 
termed  "interplanetary"  while  the  last  resides  within  the  Earth's  magnetosphere, 
and  is  termed  a  "local"  component.  The  interplanetary  particles  encompass  a  wide 
range  in  both  energy  (E  >  500  keV)  and  charge  (1  £  Z  £  92).  Some  of  these  particles 
penetrate  the  Earth's  magnetic  field  and  form  part  of  the  radiation  environment  in 
which  spacecraft  must  operate. 

For  electronic  components  located  within  a  spacecraft,  there  is  a  large  radiation 
dose  (depending  upon  the  actual  orbit)  contributed  by  protons  and  electrons  plus 
radiation  effects  due  to  the  heavy  ion  component.  This  latter  is  particularly 
important  because  of  the  ability  of  a  single  heavy  nucleus  to  affect  ("upset") 
microelectronic  circuitry  (Binder  et  al.,  1975)  particularly  low  power  memory  chips 
(Kolasinski,  et  al.,  1979).  Even  though  the  relative  abundance  of  heavy  ions  is  small, 
compared  to  protons,  the  heavy  ion  efficiency  is  high  since  the  effects  are  caused 
directly  by  the  passage  of  the  particle  through  the  device.  The  major  concern  for  the 
future  lies  in  the  decreasing  feature  size  of  the  electronic  systems  being  designed  or 
proposed,  since  as  the  feature  size  decreases,  the  sensitivity  of  these  systems  to 
disruption  by  individual  particles  increases  (see,  for  example  Iverson,  1979). 

An  important  objective  of  the  CRRES  program  is  to  investigate  the  effects  of 
the  space  radiation  environment  on  the  modern  microelectronic  components  that 
will  be  used  in  future  spacecraft.  One  of  these  effects  is  the  ability  of  single, 
intensely-ionizing  particles  to  upset  the  logical  state  of  a  single  bit  in  a  digital 
microcircuit,  causing  a  single  event  upset  (SEU).  The  MicroElectronics  Package 
(MEP)  onboard  the  CRRES  satellite  is  measuring  SEU  rates  in  a  large  sample  of 
digital  components. 

The  data  from  the  CRRES  particle  sensors  provides  simultaneous 
measurements  to  characterize  the  environment  being  sampled  by  the  engineering 
experiments.  However,  detailed  characterization  requires  significant  amounts  of 
analysis  effort  and  extended  periods  of  time,  thereby  affecting  the  ability  of  the 
engineering  experiments  to  complete  their  analysis  schedule.  An  accurate  model  of 
the  heavy  ion  radiation  environment  is  required  for  the  analysis  of  the  single-event  » 
upset  data  from  the  CRRES  Microelectronics  Package  and  to  support  the  other 
engineering  experiments  on  the  CRRES  mission.  Of  longer  range  concern  is  the 
ability  of  ground  testing  and  calculations  to  predict,  accurately,  the  effects  of  the 
space  radiation  environment  on  satellite  systems  and  instrumentation.  One  of  the 
important  goals  of  the  MEP  program  is  to  develop  such  methods  for  spacecraft 
electronic  systems.  An  essential  part  of  such  calculations  is  an  accurate  model  of  the 
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space  radiation  environment  that  is  responsible  for  the  SEU's,  particularly  the  heavy 
ion  component.  Thus,  the  predictive  model  being  developed  under  tliis  contract 
will  retain  its  usefulness  long  after  the  CRRES  mission  is  completed. 

b.  Technique 

A  computer  model  of  the  Interplanetary  Heavy  Ion  Radiation  Environment, 
calibrated  and  adjusted  with  the  actual  CRRES  flight  data  from  a  number  of  different 
particle  sensors,  which  is  updated  at  regular  intervals  throughout  the  CRRES 
mission  lifetime,  can  provide  the  data  needed  for  the  interpretation  of  the  CRRES 
Engineering  experiments,  particularly  the  Microelectronics  package,  while 
maintaining  faithfulness  to  the  actual  conditions  encountered  by  the  CRRES 
spacecraft.  The  model  under  development  is  an  analytic  description  of  the  f  \x  of 
each  ion  species  as  a  function  of  energy  and  time.  It  will  include  the  two  major 
sources  of  heavy  ion  flux  in  the  interplanetary  medium:  galactic  cosmic  rays  and 
solar  energetic  particles.  The  model  will  cover  the  energy  range  above  ~10 
MeV/nucleon  and  will  describe  all  elements  including  estimates  for  those  elements 
that  are  no:  well  measured.  Once  developed,  the  model  will  have  a  minimum 
number  of  parameters  needed  to  describe  the  overall  environment.  To  keep  the 
model  accurate  and  to  insure  that  time  variations  do  not  become  large,  periodic 
updates,  approximately  semi-annually  throughout  the  CRRES  mission  will  be 
issued.  These  updates  will  be  keyed  to  the  model  parameters  and  will  provide  the 
information  necessary  for  a  user  to  update  his  calculations  for  the  particular  period 
under  consideration.  Following  termination  of  the  CRRES  mission,  a  final  model 
release  may  be  made  along  with  the  techniques  necessary  to  utilize  and  update  the 
model  software  for  future  time  periods.  Such  a  new  predictive  model  will 
supplement  the  products  of  the  Microelectronics  test  package,  so  as  to  produce  a 
standard  method  of  estimating  SEU  rates  for  future  space  missions. 

It  might,  at  first,  seem  best  to  employ  the  measurements  from  the  particle 
sensors  on  CRRES  directly  to  determine  the  radiation  environment.  This  is  not 
practical  for  several  reasons:  (1)  the  total  flux  of  particles  is  so  low  that  accurate 
measurements  of  all  particle  species  will  take  a  long  time  to  accumulate;  (2)  a 
thorough  analysis  of  the  data  to  determine  accurate  abundances  of  each  individual 
element  will  be  time-consuming;  (3)  the  CRRES  sensors  can  only  measure  about 
half  of  the  integral  flux  during  quiet  periods  because  they  do  not  measure  all  of  the 
highest  energy  particles;  and  most  importantly,  (4)  the  scientific  research  in  cosmic 
ray  physics  over  the  past  30  years  has  established,  accurately,  relative  abundances 
and  energy  spectra  for  most  of  the  elements.  The  best  GCR  model  should  take  full 
advantage  of  this  accumulated  knowledge  and  use  contemporary  measurements 
most  efficiently  to  accurately  define  the  cosmic  ray  intensities  over  the  energy  ranges 
where  temporal  variations  are  large  and  not  easily  predicted  in  detail. 

For  these  reasons,  a  different  approach,  is  undertaken,  namely  a  description 
based  on  the  well-established  physical  principles  of  propagation  and  solar 
modulation  that  control  and  determine  the  time-  and  energy-dependent  GCR 


irensity  at  earth.  This  model  will  incorporate  obtei  vatious  obtu^ea  uvi.  tne  past 
v  :;:ades  of  space  research.  Even  when  complete,  the  CRRES  particle  sensor 
.  r  are  unlikely  to  he  a  significant  improvement  over  the  existing  GCR  database 
n  terms  of  composition,  characteristics  of  energy  spectra,  or  range  of  intensity 
venation.  Thus,  the  CRRES  particle  data  sets  will  not  provide  information  that 
would  charge  the  known  qualitative  features  of  the  interplanetary  high  energy 
heavy  ion  environment  model,  but  are  of  crucial  importance  in  obtaining 
quantitative  agreement  between  the  model  and  the  exact  conditions  during  the 
period  of  the  CRRES  mission. 

Quantitative  values  of  parameters  in  any  heavy  ion  computer  model  have  a 
significant  range  of  variation  with  time.  These  variations  are  dominated  by  solar 
phenomena  and  the  changing  levels  of  solar  activity.  While  the  solar  cycle  period 
and  phase  are  known,  active  solar  processes  are  sufficiently  random  that  short  term 
effects  cannot  be  accurately  predicted  in  advance.  CRRES  data  will  be  analyzed  on,  at 
least,  monthly  time  scales  and  will  be  subject  to  such  fluctuations.  In  addition,  for  a 
3-5  year  mission,  CRRES  should  sample  conditions  from  solar  maximum  through 
solar  minimum,  and  a  viable  model  must  be  able  to  describe  the  full  range  of  solar 
cycle  conditions.  Precise  values  for  the  parameters  of  the  model  must  depend  upon 
contemporary  measurements  of  the  heavy  ion  environment.  Particle  flux 
measurements  from  sensors  aboard  CRRES  are  thus  required  to  normalize  the 
model  as  well  as  to  describe  the  conditions  following  solar  energetic  particle  (SEP) 
events.  The  most  accurate  and  effective  way  to  merge  the  CRRES  particle  sensor 
datasets  with  the  environment  model  will  be  to  use  the  data  to  define  key 
parameters  of  the  model.  For  the  galactic  cosmic  ray  (GCR)  component  these  key 
parameters  will  describe  the  level  of  solar  modulation  and  presence  of  any 
"anomalous"  components.  For  the  solar  flare  component  these  key  parameters  will 
define  the  occurrence  of  the  flare,  its  peak  flux,  total  fluence,  composition 
characteristics,  and  the  form  of  its  energy  spectrum. 

c.  Heavy  Ion  Components 

The  model  will  provide  a  description  of  the  flux  of  each  ion  species  as  a 
function  of  energy  and  time,  focusing  on  the  two  major  sources  of  heavy  ions,  GCR 
and  SEP.  The  model  addresses  the  energy  range  above  ~10  MeV/nucleon  and  will 
describe  all  elements,  including  estimates  for  the  rare  species  that  are  not  well 
measured. 

i.  Galactic  Cosmic  Rays 

Figure  1  shows  the  relative  event  rate  for  each  of  the  elements  Z  >  2  in 
the  galactic  cosmic  rays,  with  the  points  connected  by  lines  within  each 
element  giving  the  relative  abundance  of  the  isotopes  of  that  element.  The 
point  to  note  here  is  that  the  heavy  ion  flux  is  dominated  by  the  major 
species,  C,  O,  Ne,  Mg,  Si,  Fe,  but  the  other  elements,  taken  together,  still 
constitute  a  significant  part  of  the  total  flux.  The  analysis  of  CRRES  particle 
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Figure  1.  Relative  event  rates  for  GCR  Heavy  Ions  at  solar  minimum.  The  lines 
connect  the  different  isotopes  of  the  elements. 


^  % 
Oxygen 


The  University  of  Chicago 
IMP  8 

Differential  Energy  Spectra 
1974-1976 


Magnesium 


Kinetic  Energy  (MeV/nucieon) 


Kinelic  Energy  (MeV/nucleon) 


Figure  2.  Measured  energy  spectra  for  C,  O,  Ne,  Mg  and  Si  from  IMP-8. 


sensor  data  will  focus,  initially,  on  the  most  abundant  spedes  (H,  He,  C,  O,  Ne,  Si, 
Ca,  and  Fe).  The  data  reduction  efforts  will  produce  fluxes  for  each  of  these  in 
energy  intervals  as  appropriate  and  convenient  for  each  of  the  sensors,  with 
attention  paid  to  instrument  calibrations  and  backgrounds  in  order  to  obtain 
accurate  flux  normalizations. 

Figure  2  shows  an  example  of  measured  GCR  heavy  ion  spectra  for  the 
major  species,  C,  O,  Ne,  Mg,  Si  from  below  100  MeV/ nucleon  to  about  1 
GeV/nudeon  taken  at  solar  minimum  modulation  conditions  (Garcia- 
Munoz  et  al.,  1 977).  The  spectra  peak  at  several  hundred  MeV/nudeon  and 
fall  off  on  both  the  higher  and  lower  energy  sides  due  to  the  nature  of  the 
energy  (rigidity)  spectrum  and,  at  low  energies,  the  effect  of  solar  modulation. 

The  basic  shape  of  the  individual  element  spectra  are  similar,  however,  the 
shape  of  the  spectra  do  vary  with  the  level  of  solar  modulation. 

The  cosmic  rays  diffuse  into  the  Heliosphere  against  the  outward 
flowing  solar  wind.  In  this  solar  modulation  process,  the  cosmic  rays  lose 
energy  to  the  expanding  field  (adiabatic  deceleration)  and  their  energy  spectra 
are  modified.  Thus,  the  spectra  observed  at  Earth  are  not  the  same  as  the 
spectra  outside  the  Heliosphere  in  local  interstellar  space  (LIS).  In  addition, 
the  observed  spectra  vary  with  time  as  conditions  change  within  the 
Heliosphere.  The  time  scales  for  significant  changes  in  the  intensities  of  the 
GCR  flux  are  typically  months.  Thus  a  model  of  the  galactic  cosmic  ray 
component  can  be  kept  accurate  and  up-to-date  with  only  periodic  changes  to 
its  parameters.  (Solar  modulation  is  described  in  more  detail  in  a  subsequent 
section.) 

Figure  3  shows  a  compilation  of  data  (Simpson,  1983)  for  the  energy  spectra  of 
Hydrogen,  Helium,  Carbon  and  Iron.  Above  1-10  GeV /nucleon  (depending  upon 
the  element),  the  spectra  are  power  laws.  Below  these  energies  the  spectral  shape 
varies  with  the  modulation  level.  The  CURES  sensors  measure  below  ~1 
GeV/nucleon,  and,  therefore,  the  flux  for  the  highest  energy  particles  must  be 
calculated  from  the  "known"  characteristics  of  the  spectra. 

Figure  4  shows  a  set  of  recent  measurements  (1987)  on  the  spectra  of  protons 
and  helium  nuclei  (Seo  et  al.,  1991).  The  dashed  curves  show  the  results  of 
calculations  for  different  levels  of  solar  modulation,  which  demonstrate  the  degree 
of  variability  introduced  into  the  spectra  by  changing  Heliospheric  conditions.  In  all 
of  the  calculations,  the  assumed  spectra  in  Local  Interstellar  Space  remained 
constant. 

Thus,  predicting  the  intensity  of  GCR's  at  Earth  requires  knowledge  of 

•  The  Local  Interstellar  Spectra  for  each  element 

•  The  level  of  solar  modulation  applicable  to  the  time 
period  under  consideration. 


Flux  (nr  sr  s  MeV/Nucleon 
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KINETIC  ENERGY  (MeV/nucleon)  KINETIC  ENERGY  (MeV/nucl«on) 

Figure  4:  Differential  energy  spectra  for  cosmic-ray  protons  (left)  and  helium  nuclei 
(right)  measured  in  1987  (filled  circles  for  LEAP,  open  circles  for  IMP  8) 
along  with  the  previous  balloon  measurements  (see  Seo  et  al.,  1991  for 
references  to  the  data).  The  dashed  curves  represent  the  local  interstellar 
and  modulated  spectra  with  different  amounts  of  modulation  indicated 
by  the  parameter  a,  local  interstellar  spectra  (no  modulation);  b,  $  = 

200  MV;  c,  0  *  400  MV;  d,  <j)  =  500  MV;  e,  <fr  =  600  MV;  f,  <J>  -  800  MV,  g,  <f>  = 
1000  MV. 

Determining  these  two  parameters  has  been  the  focus  of  much  of  our  analytical 
work  during  the  past  year. 

ii.  Solar  Energetic  Particles 

Solar  flare  events  are  semi-random  in  occurrence.  The  flare  activity  peaks  near 
solar  maximum  in  the  cycle,  but  there  are  events  scattered  throughout  the  full  solar 
cycle.  Flares  also  vary  widely  from  event  to  event,  and  this  requires  SEP  events  to  be 
described  individually.  Data  from  the  particle  sensors  onboard  the  CRRES  satellite 
and  from  other  sources  detect  SEP  events,  each  of  which  will  be  fit  to  a  SEP  model. 
The  results  are  then  added  to  the  galactic  cosmic  ray  component  to  obtain  a  complete 
description  of  the  heavy  ion  environment,  at  a  reasonable  minimum  level  of  data 
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reduction  and  analysis  effort.  Even  at  the  maximum  of  solar  activity,  there  will  be 
only  a  few  flares  per  month  that  will  produce  significant  fluxes  of  heavy  ions  at  the 
energies  of  interest  here.  Based  on  these  considerations,  semi-annual  updates 
should  suffice  to  review  the  solar  flare  component.  However,  for  particularly  active 
periods  (e.g.  the  late  March-June  1991  interval)  more  frequent  data  presentation  may 
be  desirable. 


There  have  been  a  number  of  efforts  over  the  years  to  model  SEP  spectra. 

These  have  resulted  in  the  use  of  power  laws  in  kinetic  energy  (Dietrich  and 
Simpson,  1978;  Cook  et  al.,  1984)  exponential  functions  in  magnetic  rigidity  (Freier 
and  Webber,  1963),  exponential  functions  in  momentum /nucleon  (Mewaldt  et  al., 
1984)  and  Bessel  function  fits  (McGuire  and  von  Rosenvinge,  1984)  to  the  spectra. 
The  exact  form  to  apply  seems  to  depend  upon  the  individual  flare,  the  energy  range 
of  the  fit  and  the  species  being  considered.  Figure  5,  from  Mewaldt  et  al.  (1984), 
shows  the  fits  obtained,  for  one  flare,  to  the  major  heavy  ion  species  from  carbon  to 
silicon.  In  this  particular  case,  the  exponential  form  appears  to  give  a  good 
representation  of  the  data. 
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Momentum  [(MeV/c)/  nucleon] 

Figure  5.  Solar  flare  energy  spectra. 


The  heavy  ions,  such  as  are  shown  on 
Figure  5,  are  often  enhanced  (relative  to 
normal  solar  system  abundances)  in  solar 
flares.  Often  the  enhancement  is  a  function 
of  charge,  Z.  For  heavy  nuclei  such  as  iron, 
potentially  the  most  important  for  the  SEU 
problem,  the  abundances  can  be  enhanced 
by  over  an  order  of  magnitude  relative  to 
oxygen. 

The  heavy  ion  enrichment  in  flares  of 
varying  sizes  has  been  the  subject  of 
numerous  studies  (McGuire  et  al.,  1985; 
Meyer,  1985a;  1985b;  Breneman  and  Stone, 
1985;  Mason  et  al.,  1986;  Guzik,  1988),  and 
these  models,  as  well  as  others,  must  be 
examined  to  select  a  representation  for  SEP 
events  that  describes  the  contributions  of 
these  events  to  the  various  heavy  ion 
abundances  and  spectra.  The  parameters  of 
such  a  representation  will  have  to  be 
determined  from  the  CRRES  data  for  each 
SEP  event.  However,  since  not  all  SEP'S 
will  be  large  enough  to  determine  all  of  the 
parameters  of  the  model,  we  will  rely  on 
published  data  for  earlier  SEP  events  to 
define  typical  values  that  can  be  adopted  for 
such  events. 
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The  elements  above  the  iron  peak,  Z  >  28,  are  potentially  more  damaging  to 
microelectronic  circuity  than  even  the  iron  peak  elements.  These  nuclei  are  found 
in  both  solar  flares  and  in  the  galactic  cosmic  rays  and  could  contribute  to  a  few 
isolated  events.  However,  the  normal  abundance  of  the  Z  >  28  elements  is  quite  low 
and  even  in  a  heavy  ion  enhanced  flare,  the  expected  numbers  of  such  nuclei  are 
only  a  few  tenths  of  a  percent  of  the  iron  flux.  A  similar  situation  applies  to  the 
GCR  (Binns  et  al.,  1981)  where  Zinc  is  only  10-3  of  iron  and  the  heavier  elements  are 
at  the  <10‘4  level. 


Solar  energetic  particle  events  are  unpredictable,  and  identification  of  SEP 
events,  especially  small  flares,  is  not  easy  with  the  ONR-604  sensor  data,  over  a 
limited  period  of  time.  It  will  be  necessary  to  use  data  from  other  sensors,  external 
to  the  CRRES  program,  to  unambiguously  identify  the  SEP  events.  The  monitoring 
provided  by  the  GOES  geostationary  satellites  and  by  the  IMP-8  interplanetary 
spacecraft  are  important  sources  of  SEP  information. 
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Figure  6:  GCR  energy  spectra  showing  the 
presence  of  anomalous  Nitrogen  and 
Oxygen  at  low  energies. 


iii.  The  Anomalous 

Component 

This  anomalous  component  is 
prominent  under  solar  minimum 
modulation  conditions,  at  energies 
below  ~50  MeV/nudeon  at  Earth, 
and  appears  to  involve  only  high 
ionization  potential  elements.  He,  N, 
O,  Ne  and  Ar.  Figure  6  (from  Garda- 
Munoz  et  al.,  1979)  shows  the  energy 
spectra  for  B,  C,  N  and  O  for  energies 
above  ~10  MeV/nucleon  compiled 
from  a  number  of  different 
experiments.  From  the  peaks  of  the 
spectra,  the  B  and  C  show  a  steady 
decrease  in  flux  with  decreasing 
energy  as  expected  from  the 
spectrum  and  solar  modulation 
models.  The  O  and  N  spedra, 
however,  show  a  flattening  and  a 
tum-up  at  energies  below  ~40 
MeV/nucleon.  This  is  due  to  the 
presence  of  the  so-called 
"anomalous"  component 
contributing  to  the  spectrum  at  low 
energies.  For  oxygen,  the  flux  at  10- 
20  MeV/nudeon  is  actually  greater 
than  the  peak  flux  at  several 
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hundred  MeV/nucleon.  The  anomalous  component  can  be  important  for  CRRES 
during  times  of  solar  minimum.  At  solar  maximum,  this  component  is  largely 
absent,  showing  that  it  is  affected  strongly  by  the  level  of  solar  modulation.  It  will  be 
necessary,  therefore,  to  model  and  include  the  solar  modulation  dependence  of  the 
anomalous  component  in  the  interplanetary  heavy  ion  model,  once  the  CRRES 
mission  enters  the  minimum  phase  of  the  current  solar  cycle. 

d.  Distribution  of  Effort 

Developing  the  model  described  above  involves  a  combination  of  modeling/ 
calculation  and  data  analysis/interpretation.  The  tasks  involved  include: 

-  Compilation  of  previous  GCR  and  SEP  data  and  reduction  of  this  dataset. 

-  Development  of  a  First-order  model  based  upon  the  accumulated  GCR  and 
SEP  dataset. 

-  Software  development  for  the  analysis  of  ONR-604  data  including  plotting 
software,  background  assessment,  normalization  programs,  correlation 
routines  and  the  like. 

-  Establish  techniques  to  measure  the  level  of  solar  modulation  during  the 
mission. 

-  Finalize  SEP  description  for  use  in  the  models. 

-  Support  launch  and  on-orbit  operations  of  CRRES  and  monitoring  of  the 
ONR-604  instrument. 

-  Establish  routine  processing  and  analysis  software  for  the  CRRES  data. 

-  Develop  ancillary  databases  for  comparison  and  use  in  the  modeling  effort. 

Thus,  the  effort  involves  a  combination  of  data  processing  and  analysis  with 
modeling  and  interpretation  of  both  new  and  existing  data.  All  aspects  must  be 
worked  simultaneously,  since  there  was  no  Air  Force  support  for  the  planned  pre¬ 
launch  modeling  activities. 

The  approximate  distribution  of  effort  among  the  collaborating  institutions  is 
shown  schematically  in  Figure  7.  The  data  flow  begins  at  Phillips  Laboratory  (AFGL) 
with  the  production  of  Agency  Tapes  and  proceeds  to  the  right  through  processing  at 
Chicago,  analysis  at  LSU,  model  calculations  at  LPARL  and  back  to  the  Working 
Groups  at  Phillips  Laboratory  in  the  form  of  model  predictions  of  LET  spectra  and 
SEU  estimates.  Significant  coordination  is  required  for  this  program,  but,  during  the 
first  year  of  effort,  this  approach  has  worked  reasonably  well. 
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3.  THE  ONR-604  EXPERIMENT 


The  particle  sensor  data  needed  to  refine  the  parameters  of  a  heavy  ion  model 
are  the  flux  and  energy  dependence  of  the  major  ion  species,  protons,  helium  nuclei 
and  heavier  elements.  The  intensities  of  the  heavier  elements  are  low,  and 
obtaining  a  reliable  estimate  of  their  flux  will  require  integrating  over  months,  in 
broad  energy  intervals.  In  this  analysis  care  must  be  taken  to  account  for  time 
dependent  changes  in  sensor/electronic  calibrations  and/or  backgrounds.  It  is 
prudent,  therefore,  to  analyze  elements  chosen  so  as  to  be  relevant  to  the  model  and 
to  give  reasonably  accurate  flux  measurements  over  an  approximately  three-six 
month  period,  e.g.  H,  He,  B,  C,  O,  Ne,  Si,  Ca  and  Fe-Ni,  analyzed  over  the  energy 
ranges  convenient  for  each  species.  The  "key”  experiment  for  this  heavy  ion 
analysis  is  the  ONR-604  instrument  built  by  the  Laboratory  for  Astrophysics  and 
Space  Research  at  The  University  of  Chicago  under  the  direction  of  John  A. 

Simpson.  The  experiment  was  constructed  with  NASA  support  from  the  CRIE 
project  (M.  Garcia-Munoz  and  J.  P.  Wefel,  Co-Investigators)  under  contract  NAS-2- 
24430  and  sponsored  for  flight  on-board  CRRES  by  the  Office  of  Naval  Research. 

a.  Instrument  Description  and  Operation 

The  ONR-604  investigation  was  designed  to  obtain  new  data  on  the  elemental 
and  isotopic  composition  of  high  energy,  heavy  ions  to  explore  processes  in  solar 
flares,  in  the  galaxy  and  within  our  Geospace  environment.  The  instrument  is 
optimized  for  studies  of  the  heaviest  ions,  up  through  the  iron  peak.  The  single  High 
Energy  Telescope  (HET),  shown  schematically  in  Figure  8  (Simpson  et  al.,  1985), 
consists  of  a  stack  of  fifteen  silicon  solid  state  detectors  surrounded  by  a  scintillation 
guard  counter  plus  external  passive  shielding  to  reduce  the  effects  of  low  energy 
background  radiation.  The  stack  contains  six  Position  Sensitive  Detectors  (PSD)  — 
upper  right  —  used  to  determine  the  particle  trajectory  (Lamport  et  al.,  1976)  followed 
by  eight  5  mm  thick  detectors  which  measure  the  particle  energy  loss  and  residual 
energy,  allowing  mass  and  charge  identification  by  the  AE-E  technique.  A  1000  micron 
thick  penetration  detector  (A)  is  located  at  the  bottom  of  the  stack. 

A  unique  feature  of  the  instrument  is  the  six  PSD’s,  for  which  ONR-604  is  the 
first  space  exposure  of  this  new  technology  developed  at  Chicago.  A  schematic  view 
of  a  PSD  is  shown  at  the  top  right  of  the  figure.  The  upper  surface  consists  of  many 
parallel  strips  of  gold  each  connected  via  a  wire  bond  to  a  resistive  divider  network. 

By  comparing  the  signal  from  the  resistive  divider  network  (P)  to  the  full  energy 
signal  measured  at  the  back  surface  (E),  the  strip  nearest  the  point  of  incidence  of  the 
particle  can  be  determined.  From  a  Fit  to  the  positions  in  the  six  PSD's,  the  trajectory 
of  the  particle  is  determined  with  an  angular  accuracy  of  better  than  1°,  as  validated 
by  accelerator  calibrations.  The  PSD's  allow  the  incident  trajectory  of  each  event  to 
be  determined,  providing  corrections  for  the  path  length  in  the  other  detectors  and 
point  of  incidence  on  the  detector  surface  as  well  as  allowing,  for  example,  the 
particle  pitch  angle  distribution  to  be  studied. 
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Figure  8.  The  ONR-604  Telescope  and  its  Position  Sensing  Detectors 


For  each  event  that  is  pulse  height  analyzed,  the  total  energy  deposited  in  each 
of  the  first  14  detectors  is  read  out.  This  allows  the  particle  range  ("ID")  to  be  used  to 
subdivide  the  full  energy  interval.  For  high  precision  charge  and  mass 
identification,  the  particle  must  come  to  rest  within  the  instrument,  which 
determines  the  energy  over  which  each  type  of  particle  can  be  studied.  The  entire 
stack  is  surrounded  by  a  scintillator  guard  counter  to  eliminate  out-of-geometry 
events  and  particles  that  exit  from  the  side  of  the  apparatus.  The  scintillator  signal 
is  not  part  of  the  trigger  logic  but  is  recorded  (along  with  A)  as  a  flag  for  each 
analyzed  event.  In  addition,  a  ground-commandable  operations  mode  is  provided 
for  the  detailed  analysis  of  proton  and  alpha  particle  flux  levels. 

Because  the  anti-coinddence  srintillator  guard  subtends  a  large  solid  angle,  it 
must  be  shielded  from  the  high  counting  rates  induced  by  electrons  trapped  in  the 
magnetosphere.  Passive  shielding  to  eliminate  <4  MeV  electrons  is  provided  by  a 
magnesium  shield  and  Tantalum  collar.  The  lower  part  of  the  scintillator  is  located 
well  within  the  instrument  package  and  is  surrounded  on  five  sides  by  electronic 
modules  which  offer  additional  shielding.  The  backward  moving  particles  are 
shielded  additionally  by  the  spacecraft  structure. 

In  order  to  protect  further  against  accumulated  radiation  damage  to  the 
detectors  Dl,  D2,  and  D3  (where  the  average  geometrical  factor  is  87  cm2-sr),  a  3  mm 
external  aluminum  window  shield  is  placed  over  the  telescope.  This  shield  lowers 
the  magnetospheric  integrated  proton  and  electron  radiation  doses  for  detectors  Dl, 
D2,  and  D3  below  the  tolerance  threshold  of  these  detectors  while  increasing,  only 
slightly,  the  energy  intervals  for  measurement  of  the  heavy  particles  which  are  the 
main  objective  of  the  experiment. 

Separating  isotopes  for  elements  as  heavy  as  nickel,  requires  precision 
measurements  of  the  energy  deposits  in  the  detectors.  For  this,  exceptionally  linear, 
highly  stable  electronic  circuits  to  analyze  the  signals  have  been  developed.  Each 
analysis  chain  is  individually  calibrated  preflight.  The  performance  of  each  analysis 
chain  is  monitored  during  flight  by  an  on-board  calibrator  which  sends  a  sequence  of 
2000  pulses  spanning  the  entire  range  of  signal  sizes.  Further,  in-flight  calibration 
and  normalization  is  provided  by  measurement  of  actual  galactic  cosmic  ray  events. 

In  addition  to  heavy  particles,  the  instrument  measures  protons  and  alpha 
particles  over  the  energy  range  25-120  MeV/nucleon.  In  quiet  regions  of  space,  a 
priority  system  insures  that  the  heavy  nuclei  are  analyzed  in  preference  to  protons 
and  alphas.  In  regions  of  high  background  radiation,  a  commandable  heavy  particle 
"normal  mode"  is  invoked.  This  raises  the  thresholds  on  the  detectors  Dl  through 
D6  to  suppress  the  high  firing  rates  induced  by  trapped  protons  and  electrons  while 
still  permitting  the  rare  heavy  ions  to  be  analyzed.  In  "normal  mode,"  the  "proton 
mode"  thresholds  of  detectors  Dl  through  D6  are  increased  by  a  factor  of  5.8  and,  as  a 
consequence,  the  instrument  becomes  practically  insensitive  to  protons. 


Figure  9.  ONR-604  energy  ranges  for  selected  isotopes. 

The  energy  ranges  for  triggering  the  different  detectors  for  different  species  in 
the  ONR-604  instrument  are  illustrated  in  Figure  9  for  selected  isotopes  from  H 
through  Ni.  The  lines  connect  the  incident  energies  for  0°  and  30°  angles  of 
incident.  Figure  9  indicates  that  H  and  He  are  studied  over  an  approximate  interval 
of  25-120  MeV/nudeon  while  Fe  is  measured  over  a  range  of  ~100-550 
MeV/nudeon.  Comparison  with  Figure  3  shows  that  these  intervals  are  just  below 
the  peak  in  the  cosmic  ray  spectra. 

The  instrument  acceptance  is  also  a  function  of  energy.  For  the  lowest  energy 
events,  those  that  stop  in  D1  -  D3,  the  geometrical  factor  is  >80  cm2-sr  while  for  the 
highest  energies  (stopping  in  K8)  the  geometrical  factor  is  only  ~3  cm2-sr.  This 
changing  acceptance  provides  one  of  the  principal  normalization  problems  in 
deriving  particle  spectra. 

The  highest  quality  data  is  obtained  for  partides  stopping  in  K2-K8.  However 
the  HET  also  measures  particles  that  (1)  stop  in  Kl,  (2)  stop  in  D2  -  D6  and  (3) 
penetrate  K8.  For  (1)  and  (2),  it  is  necessary  to  use  the  signal  from  the  E  contact  of 
the  PSDs  to  provide  energy  loss  information.  For  events  stopping  in  D2  or  D3,  the 
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path  length  uncertainty  is  given  by  sec  9  where  0  is  the  maximum  angle  through 
the  window  (~69°  for  D2)  which  gives  a  spread  of  a  factor  of  2.8.  Thus,  Dl,  D2  and  D3 
analysis  can  separate  protons  and  helium,  but  no  heavier  elements.  For  particles 
stopping  in  D4,  the  sec  0  spread  is  reduced  to  a  factor  of  1.44.  Particles  stopping  in  D5 
-  Kl,  have  trajectory  information  available,  and  elemental  separation  is  limited  only 
by  the  quality  of  the  E  signals  from  the  PSDs.  Finally,  the  charge  composition  of 
events  penetrating  the  entire  stack  can  be  studied  over  a  limited  energy  interval. 

The  ONR-604  instrument  incorporates  a  priority  system  for  the  pulse  height 
analyzed  events.  Each  particle  is  defined  by  the  instrument  logic  hardware  as  either 
a  PI,  P2  or  P3  event.  P3  is  the  instrument  trigger,  defined  by  the  logic  condition 

P3  =  GC  =  D1D2D3S  +  D1D2D3(K1  +  K2). 

Thus,  allowed  events  must  trigger  the  first  three  PSD's,  and  not  go  out  the  sides 
triggering  the  scintillator  OR  must  trigger  the  first  three  detectors  and  one  of  either 
Kl  or  K2,  independent  of  the  scintillator.  (Commands  are  available  to  allow  any  or 
all  of  detectors  Dl,  D2,  D3,  Kl,  K2,  S  or  A  to  be  disabled  in  the  logic.)  P2  events  add  to 
P3  the  requirement,  (KIM  +  K2M)A  ,  i.e.  the  particle  must  stop  in  the  stack  and 
trigger  the  medium  level  discriminator  threshold  on  detectors  Kl  or  K2.  Finally,  PI 
events  add  to  P3  the  requirement  that  one  of  the  high  discriminator  thresholds  on 
Kl  -  K8  must  be  triggered.  Thus,  considering  tne  discriminator  thresholds,  PI 
responds  approximately  to  >Ne  (stopping  or  penetrating);  P2  records  He  to  ~Ne 
(stopping),  and  P3  encompasses  everything  else  allowed  by  the  trigger  logic.  When  a 
high  priority  event  is  detected,  it  will  cancel  the  pulse  height  analysis  of  a  lower 
priority  event.  Thus  a  PI  will  replace  a  P2  or  P3  event,  while  a  P2  will  override  a  P3 
event.  (Discriminator  patterns  indicate  when  such  an  override  has  taken  place.) 

This  priority  system  ensures  that  a  maximum  number  of  PI  events  will  be  analyzed 
and  that  the  preponderance  of  H  and  He  will  not  "sv.  amp"  the  heavy  ions.  This  is 
particularly  important  for  obtaining  heavy  ion  fluxes  for  the  model,  since  for  Iron, 
for  example,  over  the  energy  range  180-535  MeV/nudeon,  corresponding  to  IDK1 
through  IDK8  in  ONR-604,  the  expected  number  of  nudei  is  only  about  a  dozen  per 
day! 


The  ONR-604  instrument  returns  both  pulse  height  and  count  rate  data,  the 
latter  involving  b^th  singles  and  coinddence  rates,  including  the  PI,  P2  and  P3 
priority  rates.  The  counting  rates  are  important  in  their  own  right,  provide  a 
diagnostic  for  instrument  performance  and  are  used  to  normalize  the  pulse  height 
readouts  to  obtain  particle  flux.  The  rates  are  accommodated  for  a  fixed  time  period 
of  4.096  seconds  and  their  readout  is  controlled  by  a  commutation  system.  Logic 
rates  are  output  continuously  while  detector  singles  counting  rates  are  sub- 
computated  so  that  singles  rates  are  available  less  frequently.  Table  1  gives  the  logic 
for  the  different  rates  as  well  as  an  indication  of  the  approximate  particle  energy 
range  sampled  by  the  rate. 
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TABLE  1:  ONR-604  Counting  Rates 


Rate  No. 

Detector 

Interval 

(Rate 

Name) 

Rate  Logic 

Energy  Interval 
of  Rate  in 

Proton  Mode 
(MeV) 

Particle 
Defining 
Energy  Interval 

1 

D1 

D1D2D3 

(S)a 

24.8  -  44.5 

Protons 

2 

D2 

D1D2  D3  D4 

(S)a 

28.3  -  48.2 

it 

3 

D4 

D13D4  D5  D6 

(S)a 

35.1  -  43.9 

it 

4 

D5 

D13D5  D6  K1 

(S)a 

36.6  -  45.5 

it 

5 

K1 

D13K1  K2  K3 

(S)a 

40.0  -  59.4 

it 

6 

K2 

D13K1K2  K3 

(S)a 

52.0  -  70.6 

it 

7 

K3 

D13K1K3  K4 

(S)a 

62.1  -  79.9 

ti 

8 

K4 

D13K1K4  K5 

(S)a 

71.2  -  88.5 

i* 

9 

K5 

D13K1K5  K6 

(S)a 

79.3  -  96.5 

it 

10 

K6 

D13K1K6  K7 

(S)a 

87.0  - 103.5 

it 

11 

K7K8 

D13K1K7A 

(S)a 

94.0-116.9 

it 

12 

A 

D13K1K8A 

(S)a 

107.5-275 

♦i 

13 

P3 

D13S+D13(K1+K2) 

24.8  -  275 

it 

14 

P2 

P3  (K1M+K2M)A 

42- 105 

He4 

15 

PI 

P3  (K1H+K2H+K3H+ 
K4H+K5H+K6H+K8H) 

127  -  352  (typical) 

Mg24 

16 

Single's  Rate  Sub-corn 

Notes:  M  =  Medium  Threshold;  H  =  High  Threshold;  (Low  Threshold  is  indicated 
by  detector  symbol) 

S  Scintillation  guard  counter 

(S)a  S  term  included  every  alternate  readout. 

P3  Pulse  Height  Analysis  Gate  Control  (Priority  No.  3  or  low  priority). 

P2  Priority  No.  2  or  Medium  Priority 
PI  Priority  No.  1  or  High  Priority 


18 


Figure  10  shows  a  sample  plot  of  the  rates  for  orbit  649,  first  part  of  day  108  of 
1991.  Figure  10a  gives  the  three  priority  rates  and  the  command  state  of  the 
instrument.  ONR-604  went  into  proton  mode  (PM)  about  0615  and  remained  there 
until  -1115,  corresponding  to  the  apogee  portion  of  the  orbit.  The  P3  rate,  for 
example,  jumped  about  an  order  of  magnitude  when  PM  mode  was  on.  Note  the 
PI,  P2  and  P3  spike  around  0900.  This  corresponds  to  a  calibration  period  as 
indicated  by  the  CAL  (C)  command,  and  must  be  removed  from  rate  averages. 

There  was  also  a  data  gap  during  this  orbit  from  -1145  - 1350.  Just  after  the  gap  the 
P2  rate  shows  a  tremendous  "spike"  for  a  single  readout.  These  occur  periodically 
and  must  also  be  eliminated  to  obtain  clean  data.  Finally,  one  should  note  the  PI 
rate  in  which  the  individual  particles  can  (almost)  be  counted,  indicating  the  low 
intensity  of  the  heavy  ions. 

Figure  10b  displays  the  singles  rates  from  the  different  detectors  as  -1  minute 
averages,  and  Figures  lOc-e  show  the  coincidence  rates  listed  in  Table  1.  In  each  plot 
of  Figure  lOc-e,  the  rate  logic  with  theS  term  is  plotted  as  the  solid  line,  while  the 
rate  without  this  term  is  given  as  the  dashed  line.  Rate  D1  and  D2  are  shown  as  33 
second  averages  while  the  remainder  are  displayed  as  3  minute  averages.  In  all 
cases,  the  peak  between  0400  and  0500  corresponds  to  low-L  values  (L  =  1-3)  and 
shows  the  effect  of  magnetospheric  particles  on  the  ONR-604  counting  rates.  Note 
also  that  there  is  significant  structure  in  the  rates  during  a  perigee  pass  (0400-0500) 
and  into  the  "apogee"  portion  (0630,  L>5)  of  the  orbit.  Understanding  this  structure 
is  one  of  the  challenges  for  the  ONR-604  data  analysis. 


b.  Pre-launch  Calibration 

The  ONR-604  instrument  was  calibrated  in  July,  1977,  April,  1978,  July,  1979, 
and  again  in  February,  1989  with  beams  of  He,  C,  Ne,  Ar,  and  Fe  particles  (with 
energies  such  that  the  particles  stopped  in  the  telescope)  at  the  Bevalac  heavy  ion 
accelerator  at  the  Lawrence  Berkeley  Laboratory.  These  extensive  calibrations 
demonstrated  that  ONR-604  can  resolve  isotopes  of  H  to  Ni,  over  the  relevant 
energy  ranges.  In  addition,  these  runs  provided  a  means  to  check  instrument 
operation,  priorities,  rates  and  thresholds.  Finally,  the  use  of  a  radioactive  source 
during  the  Feb.  89  calibration  provided  a  simulation  of  expected  instrument 
performance  in  a  high  radiation  background  environment. 

For  the  acceleration  work,  special  interface  electronics  were  constructed.  This 
interface  eliminated  the  slow  synchronous  readout  of  the  instrument  and  replaced 
it  with  a  fast  asynchronous  readout  which  allowed  the  accumulation  of  significant 
statistical  samples  of  events  from  the  accelerator  beam.  For  final  confirmation  of 
performance,  the  CRRES  S/C  interface  was  re-installed  and  data  were  recorded  in 
spacecraft  mode.  Significantly  more  heavy  ion  events  were  recorded  during  these 
accelerator  runs  than  will  be  seen  during  the  entire  CRRES  mission! 


> 
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Phoenix -2  Orbit  649  Louisiana  State  University 

(CRRES  ONR-604)  Pri.  Rates  +  Cmds  Begins:  D108,  1991  (18  Apr)  Plotted:  12  Dec  91 
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Figure  10a.  ONR-604  Priority  Counting  Rates  and  Command  State. 


Phoenix -2  Orbit  649  Louisiana  State  University 

(CRRES  ONR -604)  Singles  Rates  Begins:  0108.  1991  (18  Apr)  Plotted:  12  Dec  91 
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1991 


’hoenix-2  Orbit  649  Louisiano  State  University 

CRRES  ONR -604)  KS-A  Rates  Begins:  0108,  1991  (18  Apr)  Plotted:  12  Dec  91 
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For  the  accelerator  calibrations,  the  instrument  was  mounted  in  a  gimballed 
cradle  that  allowed  it  to  be  rotated  in  two  axes.  This  permitted  the  beam  to  enter  the 
telescope  at  a  variety  of  zenith  and  azimuth  angles.  Such  angular  variations 
allowed  the  performance  of  the  PSD  trajectory  system  to  be  checked  over  the  full 
range  of  allowed  particle  incidence  angles  and  permitted  the  study  of  "edge  effects" 
for  particles  incident  at  the  maximum  possible  angles. 

The  accelerator  beams  were  first  focussed  at  the  instrument.  Then,  in  front  of 
the  instrument  was  placed  a  polyethylene  target  (to  fragment  the  beam  into  lower 
charge  isotopes)  and  a  rotatable  aluminum  wedge.  The  rotation  of  the  wedge 
reduced  the  energy  of  the  beam  to  provide  events  stopping  throughout  the  detector 
stack.  The  beam  was  then  defocussed,  which  expanded  it  from  about  1  in2  to  ~  4  in2, 
to  cover  the  full  size  of  the  front  three  detectors  (see  Figure  8).  Small  adjustments  in 
beam  energy  were  made  by  a  remote  controlled  absorber  system,  in  order  to  study 
events  stopping  in  the  PSD's  or  near  the  top  and  bottom  surfaces  of  one  of  the  K 
detectors. 


For  some  runs  it  was  necessary  to  generate  the  ion  of  interest,  e.g.  helium.  This 
was  done  by  putting  a  thick  target  in  the  beam  upstream  of  the  bending  and  focusing 
magnets.  The  isotope  of  interest  was  then  focussed  by  the  magnets  onto  the 
telescope  (of  course,  all  isotopes  with  the  same  rigidity  are  focussed  as  well,  so  these 
generated  beams  are  not  as  clean).  This  allowed  a  specific  charge  range  to  be 
investigated  in  detail. 


Figure  11:  Raw  Fe  and  Ne  calibration  data  from  the  Bevalac  showing  the  intensity 
variations  with  energy. 


The  energy  coverage  obtained  with  the  wedge  system  was  not  uniform  as 
shown  in  Figure  11  for  Iron  plus  Neon  data.  The  large  black  blobs  correspond  to 
concentrations  of  beam.  Thus,  in  the  analysis  it  was  necessary  to  remove  events  to 
try  to  obtain  as  uniform  a  track  coverage  as  possible,  within  the  constraint  of  still 
maintaining  ample  statistics. 

Figure  12  displays  samples  of  the  data  from  three  separate  "beams,"  Iron  (top). 
Neon  (center)  and  Helium  (bottom).  In  all  cases  the  plot  shows  the  energy  deposited 
in  D1  +  D2  vs  the  sum  of  the  energy  deposits  in  D3  through  K8.  For  the  Iron  beam 
only  priority  1  events  are  shown.  By  combining  data  from  several  beams  it  has  been 
possible  to  populate  the  element  "tracks"  over  the  full  matrix.  Note  the  tracks  that 
are  almost  vertical,  labeled  "R"  on  the  Neon  plot.  These  correspond  to  events  that 
passed  through  D1  D2  D3  but  missed  the  edge  of  D4  or  events  that  passed  through 
D1-D6  but  missed  the  edge  of  Kl.  Recalling  the  geometry  of  Figure  8,  the  detector 
radii  decrease  between  Dl-3,  D4-6  and  Kl-8,  so  such  events  are  allowed.  Some  will 
pass  into  the  scintillator  guard  counter  (S),  but  others  will  stop  in  the  structural 
materials  before  S.  This  is  one  source  of  "background"  that  must  be  eliminated 
from  the  flight  data. 

There  is  also  a  black  horizontal  track  on  the  Neon  plot  which  corresponds  to 
events  leaving  the  edge  of  the  stopping  detector  and  not  penetrating  to  S.  Such 
events  can  be  eliminated  by  using  the  instrument  trajectory  system  to  require  the 
particle  to  be  incident  upon  the  central  2  cm  of  the  last  detector  triggered.  The  result 
of  such  a  "clean-up"  is  shown  in  the  top  panel  of  Figure  13,  and  has  already  been 
applied  to  the  Fe  and  He  parts  of  Figure  12.  Note  also  that  one  can  begin  to  see 
separated  isotopes  for  H  and  He  and  the  fragments  of  Ne. 

Figure  13  shows  the  high  resolution  data,  events  that  stop  in  detectors  K2-K8,  as 
the  energy  deposited  in  Kl  versus  the  sum  of  energies  in  K2-K8,  for  the  Neon 
calibration.  Note  that  distinct  isotope  tracks  for  the  Neon  fragments  are  visible  The 
lower  portion  of  Figure  13  shows  the  data  restricted  to  P2  and  P3  events,  i.e.  PI  not 
triggered.  The  high  discriminator  thresholds  on  the  K  detectors  (used  to  determine 
the  P1-P2  division)  cross  the  Neon  energy  deposit  curve  dividing  this  element 
between  priority  1  and  priority  2  as  shown.  Such  a  calibration  allowed  the  KH 
discriminator  thresholds  to  be  checked  and  compared  to  the  design  values  and 
calculated  track  structure. 

Figure  14  shows  a  combination  of  the  Ne  and  Fe  datasets  with  the  Iron  track 
eliminated  so  that  the  fragment  tracks  can  be  observed.  Over  portions  of  the  energy 
interval,  separate  isotope  tracks  can  be  discerned  up  into  the  sub-iron  region.  This 
demonstrates  the  intrinsic  ability  of  the  instrument  to  resolve  heavy  isotopes.  In 
the  lower  portion  of  the  figure,  the  data  are  restricted  to  only  PI  events.  The 
portions  of  the  Neon  track  in  PI  are  apparent.  Closer  inspection  shows  that  F  and 
Na  are  also  divided  between  PI  and  P2  as  predicted  by  the  calculations. 
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Figure  12:  Accelerator  calibration  data  from  Iron  (top).  Neon  (center)  and  Helium 
(bottom)  runs. 


K 1  (MEV)  K 1  (MEV) 


Figure  13.  Bevalac  Neon  data  showing  all  events  (top)  and  priority  2+3  events 
only  (bottom).  A  stopping  radius  cut  has  been  applied  to  the  data. 
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Looking  at  the  Helium  events  allows  the  P2-P3  transition  to  be  investigated  as 
well  as  the  instrument  command  states.  Figure  15  shows  plots  similar  to  Figure  1? 
Uicttom).  In  the  top  panel  the  instrument  is  in  its  "normal"  mode  command  state 
(i.e.  high  D  thresholds  to  eliminate  protons  and  electrons)  to  be  compared  to  the 
'"proton"  mode  state  which  was  shown  in  Figure  12.  Note  that  all  of  the  hydrogen 
isotopes  are  eliminated  by  the  normal  (=heavy)  mode  command.  The  lower  part  of 
Figure  15  shows  data  restricted  to  PI  and  P2  events  (i.e.  no  P3  triggers).  This 
indicates  that  the  majority  of  the  helium  track  lies  in  P2,  but  the  lowest  and  highest 
energy  helium  events  are  found  in  the  P3  priority  state.  Thus,  reconstructing  the 
Helium  (Neon)  energy  spectrum  will  require  normalizing  the  P3  (PI)  priority 
readout  to  the  P2  readout  efficiency. 


As  part  of  the  1989  calibration,  we  investigated  the  operation  of  the  instrument 
in  a  high  radiation  background  provided  by  a  specially  purchased  144Ce  -  144Pr 
radioactive  source  of  strength  100  milli-curies.  This  isotope,  half-life  of  284  days, 
decays  by  beta  emission  to  excited  states  of  144Pr  which  transition  quickly  to  the 
ground  state  emitting  a  gamma  ray  or  themselves  decay  by  beta  emission.  The 

maximum  beta  energy  is  just  under  3  MeV  with  a  2.2  MeV  y-ray  as  the  highest 
energy  photon.  The  source  was  prepared  on  a  type-A  mount  and  inserted  into  a 
machined  aluminum  holder  equipped  with  removeable  brass  shields,  each  of  which 
had  a  different  size  hole  machined  into  it.  This  allowed  the  intensity  of  the  source 
to  be  varied  with  a  minimum  amount  of  handling  and  no  direct  exposure. 

Data  were  taken  with  the  source  located  at  different  distances  from  the 
instrument,  using  the  full  flight  electronics.  Individual  electrons  should  not  be  able 
to  trigger  a  P3  event,  but  2  or  3-fold  electron  pile-up  could.  The  observation  of  pulse 
height  analyzed  events  produced  with  only  the  source  near  the  instrument 
confirmed  this  prediction.  Similar  "events"  will  probably  be  observed  in  the  inner 
magnetosphere  or  in  the  slot  region. 

A  second  question  was  what  effect  the  presence  of  background  radiation  would 
have  on  heavy  ions  observed  at  the  same  time.  To  answer  this  question,  the 
instrument  was  operated  with  the  Iron  beam  on  and  with  the  radioactive  source 
mounted  at  an  angle  exposing  the  front  portion  of  the  instrument.  Figure  16  shows 
the  result.  Plotted  is  the  energy  deposit  in  D1-D6  versus  the  energy  deposit  in  K1-K8 
for  an  iron  beam  run  without  the  source  (top)  and  with  the  source  (bottom).  The 
bottom  plot  shows  a  second  iron  track,  below  the  beam  track,  at  an  intensity  -1/3  of 
the  normal  track.  The  lowered  track  results  from  the  D  detector  amplifiers  being 
below  baseline  when  the  iron  event  occurred,  due  to  the  high  firing  rate  of  the  front 
detectors  from  the  source  radiation.  That  this  effect  is  limited  to  the  PSD's  is 
demonstrated  in  Figure  17  where  the  K1  vs  K2-K8  matrix  is  shown  for  the  same 
data.  No  second  track  appears  h.  the  lower  (source  "on")  plot  in  Figure  17.  This  is  as 
expected  since  the  energies  of  the  source  radiation  are  insufficient  to  cause  rapid 
firing  of  the  K  detectors.  Figure  18  shows  the  charge  histograms  in  the  two  cases 
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Figure  15.  D1  +  D2  versus  D3-K8  for  Helium  events  showing  the  effects  of  the 
normal  mode  command  (top)  and  the  P2-P3  division  (bottom). 
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Figure  17.  K1  versus  K2-K8  matrix  for  Iron  beam  events  with  NO  source  (top) 
and  WITH  source  (bottom). 
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shown  in  Figure  17,  and  there  appears  to  be  no  discernible  degradation  in  the  charge 
resolution  from  the  presence  of  the  radioactive  source.  In  addition,  there  is  no 
discernible  shift  in  the  position  of  the  Iron  track  between  the  two  panels  of  Figure 
17. 


The  majority  of  the  source  effect  is  confined  to  the  first  detector  Dl.  This  is 
illustrated  in  Figure  19  which  shows  plots  of  the  energy  deposited  in  Dl  (top),  D2 
(middle)  and  D3  (bottom)  versus  K1-K8.  Only  in  the  top  plot  is  the  "second  track" 
noticeable. 

We  expect  such  effects  to  occur  during  the  CRRES  mission,  particularly  near 
the  outer  belt.  The  flux  and  energy  distribution  of  the  background  radiation  in  space 
will  be  different  from  that  obtained  in  this  source  calibration,  however  the  effects 
should  be  predictable.  We  expect  that  the  K1-K8  events  will  be  largely  unaffected, 
but  extracting  the  lower  energies  from  the  ONR-604  data  will  require  very  careful 
analysis. 


C.  Launch  and  On-Orbit  Checkout 

Much  of  the  initial  several  months  of  effort  under  this  contract  was  devoted  to 
pre-launch  preparations,  rehearsals  and  training  at  CSTC,  launch  operations  and 
instrument  turn-on  and  initial  operations  planning  for  the  experiment.  The 
important  questions  for  ONR-604  were  the  instrument  calibration,  the  performance 
in  various  parts  of  the  orbit  in  different  command  states  and  the  development  of  a 
plan  for  routine  operations  of  ONR-604  on  CRRES.  In-flight  calibration  runs  were 
scheduled  over  various  parts  of  the  CRRES  orbit  as  well  as  while  the  instrument 
was  in  the  different  command  states,  i.e.  proton  mode  and  normal  mode.  These 
data  were  analyzed  on  a  micro-computer  system  at  CSTC  before  pronouncing  the 
experiment  "healthy"  and  implementing  the  routine  operations  plan. 

Figure  20  shows  the  results  of  an  In-Flight  Calibration  (IFC)  run  analyzed  at 
CSTC.  The  instrument  was  in  normal  mode  and  the  run  was  taken  during  a  pass 
through  the  outer  belt  region.  Plotted  are  the  deviations  of  the  IFC  points  from  a 
best  fit  straight  line.  Of  significance  is  the  spread  of  the  points  which  appears  "bad" 
for  PI,  El,  P2  and  E2,  the  amplifiers  on  the  top  two  detectors  in  the  telescope.  For  the 
other  amplifier  chains  the  deviations  appear  normal.  When  analyzed  completely, 
the  gains  deduced  from  this  run  were  all  <0.15%  different  compared  to  IFC  run 
taken  on  the  ground  well  before  launch.  Thus,  it  appears  that  the  instrument 
maintained  its  absolute  calibration  through  launch  operations.  Further,  data  taken 
in  high  background  environments  such  as  the  outer  belt  may  be  subject  to 
variations,  at  least  in  the  top  two  PDS's. 

Figure  21  shows  the  P3  counting  rate  as  a  function  of  time  for  a  full  orbit  of  data 
analyzed  at  CSTC.  The  instrument  was  in  proton  mode  at  the  time.  Note  that  P3 
responds  strongly  to  both  the  inner  and  outer  belts  but  is  relatively  quiet  under  the 
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Figure  19.  Energy  deposits  in  D1  (top),  D2  (middle)  and  D3  (bottom)  versus  K1-K8 
for  iron  events  taken  WITH  the  radioactive  source. 
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Figure  21.  One  of  the  first  full  orbits  of  data  from  ONR-604.  Plotted  is  the  P3 
counting  rate  versus  the  Vehicle  Time  Code  Word  (VTCW)  for  the 
Apogee  and  Perigee  portions  of  the  orbit  (GSE  Run  515). 
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Figure  22.  P3  and  P2  rates  versus  VTCW  for  a  full  orbit  of  ONR-604  data  in  proton 
mode  (GSE  Run  521). 
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at  perigee  and  during  this  apogee  pass.  Figure  22  shows  a  comparison  oi  the  P3 
and  P2  rates  for  a  later  pass  with  the  instrument  still  in  proton  mode.  Two  features 
« marked  A  and  B)  appear  in  the  slot  region  and  appear  to  have  longitudinal 
symmetry.  The  surprise  here  is  the  lower  plot  in  which  P2  (nominally  Z  >  2)  shows 
a  strong  response  in  the  inner  belts  but  is  insensitive  to  the  outer  belt  particle 
population.  The  inner  belt  response  may  be  to  proton  interactions  or  to  event  pile- 
up  (such  as  observed  during  the  accelerator  calibration  runs),  but  the  exact  nature  of 
the  response  must  be  investigated  further. 

The  data  shown  on  Figure  21  and  22  indicate  that  (1)  the  P3  counting  rate  (and 
the  Dl,  D2,  D3  amplifiers)  reach  exceptionally  high  values  in  the  inner 
magnetosphere,  and  (2)  the  P2  pulse  height  analysis  will  work  only  outside  of  the 
inner  zone.  Thus,  it  was  prudent  to  protect  the  top  detector  systems  by  commanding 
the  instrument  into  normal  (=heavy)  mode.  However,  this  eliminates  the 
measurements  of  protons  at  apogee.  Therefore,  an  operational  plan  was  developed 
and  implemented  at  CSTC  whereby  on  every  orbit  the  instrument  is  commanded 
into  proton  mode  at  an  altitude  of  -26,600  km  outbounds  and  is  returned  to  normal 
mode  at  the  same  altitude  on  the  inbound  leg  of  the  orbit.  This  is  the  operational 
sequence  shown  in  Figure  10a  and  has  been  followed  for  essentially  all  of  the 
mission  thus  far. 

Figure  23  shows  an  expanded  view  of  the  outer  belt  region  as  seen  in  the  D2 
and  the  P3  counting  rates.  Note  the  strong  spin  modulation  in  the  rates  (spin  rate 
-2.6  rpm  at  this  time).  Demodulating  the  rates,  in  order  to  look  at  data  on  short 
lime  scales  (-4  seconds),  will  be  necessary  to  understand  the  instrument  response 
within  the  radiation  belts. 


d.  Post-launch  Performance 

The  ONR-604  experiment  is  monitored  continuously  through  periodic  IFC 
runs,  through  the  singles  and  coincidence  counting  rates  and  through  the  pulse 
height  analyzed  events.  With  one  exception  (discussed  below)  the  instrument 
continues  to  function  nominally.  The  engineering  parameters  are  also  watched  to 
ensure  a  constant  temperature,  and  there  has  been  little  or  no  short  term  variation 
and  an  acceptable  long  term  drift. 

ONR-604  has  suffered  one  minor  degradation.  It  was  noticed  that  in  late 
October  the  PI  amplifier  failed  to  produce  a  full  signal.  The  maximum  output  was  a 
factor  of  10  smaller  than  expected.  The  other  amplifier  on  that  detector  (El) 
continues  to  respond  normally  to  both  IFC  runs  and  to  particle  events.  Thus,  we 
conclude  that  the  problem  is  not  in  the  detector  itself,  but  is  probably  a  component 
failure  in  the  pre-amplifier  circuit.  In  any  case,  the  position  information  from  Dl  is 
degraded.  Position  is  given  by  the  ratio  PI /El  so  that  a  PI  signal  down  by  a  factor  of 
10  will  give  a  position  always  on  one  side  of  the  detector.  In  fitting  a  trajectory  to  the 
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Figure  23.  Expanded  plot  of  the  P3  and  D2  counting  rates  for  an  outer  belt  pass 

approaching  Apgoee  (GSE  Run  508)  showing  the  spin  modulation  in  the  rates 
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i  cO  positions,  this  aberrant  point  will  affect  the  fit  and  will  snov.  up  o  a  iaige  at 
residual,"  a  parameter  that  is  calculated  for  each  trajectory.  However,  there  are  three 
•SD's  CD1.  D2,  D3),  for  redundancy,  in  the  top  layer,  so  that  the  loss  or  position 
information  from  one  of  the  three  is  not  a  significant  degradation  to  the 
instrument's  capability  for  reconstructing  the  particle  trajectory. 

The  difficulty  here  is  that  the  processing  software  had  been  constructed  for  six 
PSD  signals.  A  new  algorithm  had  to  be  implemented  that  used  only  five  signals. 
This  then  had  to  be  de-bugged  and  tested  and  incorporated  into  the  processing 
software.  This  halted  all  processing  until  the  algorithm  was  available,  and  required 
reprocessing  of  some  of  the  affected  orbits.  However,  when  compared  to  the 
trajectories  derived  from  six  PSD's,  the  five  PSD  results  are  essentially  identical.  It  is 
difficult  to  see  any  loss  of  precision  in  the  trajectories.  PI  is  still  being  monitored 
but,  so  far,  has  not  returned  to  good  health. 

Figure  24  shows  a  comparison  of  the  PI  events  that  have  been  analyzed  for  the 
first  ~9  months  of  the  mission  with  the  accelerator  calibration  data.  Note  that  the 
segmented  Neon  track  appears  in  both  plots.  Mg  and  Si  are  reasonably  distinct  and 
the  iron  track  appears  in  the  proper  location.  The  flight  statistics  are  still  too  meager 
to  unambiguously  identify  many  of  the  other  elements,  but  preliminary  charges  can 
be  assigned  based  upon  calculated  positions  of  the  tracks.  From  this  comparison,  we 
conclude  that  the  amplifiers  have  remained  in  calibration  and  that  the  detectors 
continue  to  function  properly. 


4.  DATA  SYSTEM  DEVELOPMENT 

In  order  to  accomplish  the  tasks  previously  discussed,  the  ONR-604  raw  data 
must  be  processed  to  (i)  convert  the  data  to  "physics"  units  (e.g.  particle  trajectory 
angles  in  degrees,  energy  deposits  in  MeV),  (ii)  obtain  auxiliary  information 
necessary  for  maintaining  the  health  and  calibration  of  the  instrument,  and  (iii) 
prepare  the  data  for  subsequent  analysis.  This  processed  data  must  then  be  further 
reduced  to  selected  rates  and  elemental  fluxes  for  particular  geomagnetic 
environments  and  times  before  information  on,  for  example,  the  level  of  solar 
modulation,  active/quiet  times  and  GCR/SEP  composition  and  spectra  can  be 
obtained.  Much  of  this  years  effort  at  LSU  has  been  associated  with  developing  such 
a  data  system,  including  the  correct  procedures  and  software  for  data  processing  and 
reduction.  To  date,  most  of  the  fundamental  processing  activities  have  been 
established,  but  additional  effort  is  still  required  to  implement  the  higher  level  data 
reduction.  In  addition,  the  system  must  include  ancillary  data  that  will  be  employed 
for  model  analysis. 
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a.  Reduction  and  Processing 


ihe  data  flow  through  the  current  lSu  processing  and  reduction  software  is 
shown  in  f  igure  2d.  Data  from  ONR-604  is  received  by  LSU  from  the  University  of 
Chicago  (UC)  in  the  form  of  a  processed  C2ET  data  tape.  These  tapes  contain  the  raw 
and  processed  PHA  events,  rate  records,  engineering/housekeeping  data  and 
associated  magnetometer  and  ephemeris  information  for  5  consecutive  orbits  with  a 
file  for  each  orbit.  Startup  delays  in  generating  Agency  Tapes  for  UC  and  final 
debugging  of  the  UC  processing  software  delayed  the  beginning  of  routine 
production  of  C2ET  tapes  until  the  beginning  of  1991.  Currently,  however,  LSU  is 
now  receiving  10-15  tapes  from  UC  every  3  to  4  weeks.  As  of  7  October  1991  LSU  has 

received  168  C2ET  tapes  from  UC  containing  data  up  to  orbit  840  (=15  July  1991)  and 
haj>  completed  initial  processing  up  to  orbit  725  (=22  May  1991).  Initial  on-orbit 
spacecraft  and  instrument  checkout  occurred  during  the  first  60-70  orbits  of  these 
data  and  correct  attitude  information  is  not  available.  Thus,  the  first  GCR  Model 
Report  will  include  6  months  of  data  beginning  ~1  September  1990. 

The  initial  step  in  the  LSU  processing  is  to  use  the  COPY_FILES  program  to 
generate  a  duplicate  of  the  C2ET  tape  in  a  format  compatible  with  the  LSU 
VAX/VMS  processing  equipment.  The  program  LISTJFILES  is  then  run  on  the 
LSU  C2ET  to  validate  the  duplicate  and  generate  an  "overview7'  listing  of  the  tape 
contents.  Once  the  LSU  C2ET  tapes  are  verified  the  original  tapes  are  returned  to  UC 
for  reuse.  Backup  versions  of  the  C2ET  are  archived  at  UC  and  only  one  copy  of 
these  tapes  are  kept  at  LSU. 

The  next  major  operation  is  to  use  the  program  RATE_GEN  to  generate  a  Rate 
Parameter  Tape.  This  tape  is  formatted  specifically  as  input  to  the  PIMPLOT 
production  plotting  program  and  contains  data  on  the  singles,  coincidence  and  PHA 
priority  rates  as  well  as  some  ephemeris  and  housekeeping  data.  Time  history  plots 
(see  Figure  lOa-e),  of  these  data  are  generated  for  each  orbit  on  a  pen  plotter,  are 
archived  and  are  used  to  monitor  instrument  health  and  calibration,  to  monitor 
geomagnetic  particle  activity,  and  for  the  initial  analysis  of  Solar  Energetic  Particle 
events. 

The  detailed  analysis  of  the  ONR-604  data  usually  requires  only  specific 
portions  of  the  full  dataset  contained  on  the  C2ET  tape.  Thus,  as  indicated  in  Figure 
25,  the  program  SPLIT_FILES  is  used  to  generate  several  data  subsets  which  are  then 
archived  for  later  use.  Processing  and  reduction  software  for  the  subsets  is  still 
under  development  and  we  have  initially  concentrated  on  interpreting  the 
instrument  rates  as  these  are  the  key  to  reliably  selecting  particle  active/ quiet  times 
and  to  determining  the  absolute  particle  flux.  As  a  result  several  programs  now 
exist  which  process  the  Rates  dataset  into  higher  order  information.  The  program 
SC_Z  provides  information  on  the  precise  position  and  orientation  of  the 
instrument  with  respect  to  the  geomagnetic  field.  This  is  needed  for  determining 
geomagnetic  cutoffs  and  for  tracing  particle  trajectories  through  the  field  in  order  to 
determine  if  the  event  is  a  trapped  particle  or  has  penetrated  from  interplanetary 
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Figure  25.  Data  Flow  through  the  current  LSU  Processing  System. 
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space.  The  program  LISTJffMES  identifies  start  and  stop  times  for  daca  gaps,  orbi.s, 
'FC  runs,  gain  changes  and  other  special  events.  This  information  is  needed  for 
correctly  determining  and  evaluating  rate  averages.  Using  SPLI  i_FiLfcS  the  Rate:, 
data  is  further  selected  for  that  portion  of  the  spacecraft  orbit  whicn  is  nearest 
apogee,  at  the  highest  L  shell  values  (generally  greater  than  4)  and  also  when  the 
instrument  gain  changes  to  "Proton  Mode".  Proton  Mode  allows  ONR-604  to  collect 
proton  events  which  are  normally  excluded  during  passes  through  the  intense 
radiation  belts.  Averages  of  these  rates  for  5  minute,  1/2  hour  and  orbit  time 
periods  are  generated  by  the  RATE_AVERAGE  program  and  time  history  plots  are 
generated  once  again  by  PIMPLOT.  These  "Proton  Mode  Rates"  are  then  used  to 
examine  the  particle  environment  which  is  least  affected  by  geomagnetic  effects  and 
which  should  be  dominated  by  interplanetary  GCR  and  SEP  particles. 

Currently,  software  to  process  the  Event  data  subsets  is  under  development. 
Two  such  subsets  are  split  from  the  C2ET  tapes  according  to  their  assigned  PHA 
analysis  priority;  Priority  1  is  assigned  to  events  which  appears  to  be  of  charge  ~1 0  or 
above  and  Priority  2  is  assigned  for  charge  ~2  and  above.  PHA  events  assigned  the 
remaining,  lowest  priority  (3)  will  also  be  split  from  the  C2ET  tapes  once  methods 
for  handling  the  electron  and  proton  pile-up  background  evident  in  these  low 
charge  data  are  developed.  Current  effort,  however,  is  focused  upon  obtaining  a 
measure  of  the  Solar  Modulation  Level  as  a  function  of  time  from  the  Priority  1  and 
Priority  2  events.  This  involves  converting  the  PHA  counts  into  a  "rate”  for  use  in 
determining  instrument  live  time  and  event  collection  efficiency;  further 
subdividing  the  data  into  quiet  and  active  times  to  separate  the  data  associated  wdth 
the  GCR  from  that  associated  with  SEP;  determining  the  event  charge,  energy  and 
trajectory  just  outside  of  the  spacecraft;  isolating  geomagnetically  trapped  and 
penetrating  particles  using  trajectory  tracing  or  cutoff  calculation;  and  applying 
appropriate  instrument  geometry  factors,  collection  efficiencies  and  live  times  to 
determine  absolute  elemental  fluxes.  The  Solar  Modulation  Level  is  then  obtained 
by  adjusting  this  parameter  in  the  GCR  Model  until  a  best  fit  is  obtained  with  the 
elemental  fluxes.  After  this  effort  is  completed  we  will  then  move  on  to 
establishing  the  processing  methods  for  modeling  SEP  events.  At  this  point  SEP 
active  times  should  already  be  isolated  in  separate  datasets  and  the  raw  data  should 
be  converted  in  flux  measurements.  Thus,  the  remaining  effort  here  will  be  to 
model  the  time  profile  and  to  determine  the  relative  enhancement  of  heavy  ions 
for  each  SEP  event. 

Analysis  of  the  ONR-604  particle  environment  data  is  neither  straight  forward 
nor  simple  due  to  the  complexity  of  both  the  instrument  and  the  environment.  It  is 
likely  that  with  improved  understanding  of  both,  the  data  processing  and  reduction 
requirements  for  ONR-604  will  be  altered.  This  can  lead  to  unanticipated 
developments  or  simplifications  to  the  processing  scheme  outlined  in  Figure  25. 


46 


b.  Ancillary  Databases 


While  the  ONR-604  data  is  the  major  component  in  the  LSU  analysis  effort,  it 
can  not  provide  all  of  the  information  required  for  a  complete  model  of  the  particle 
environment.  Thus,  LSU  has  begun  to  establish  a  database  of  auxiliary  information 
which  will  be  used  to  complement  the  ONR-604  analysis.  Components  of  this 
database  which  are  now  established  and  maintained  at  LSU,  are  shown  in  Figure  26. 

The  auxiliary  data  include  IMP-8  proton  and  helium  fluxes  and  quiet  time 
selections  in  various  energy  bins  and  the  Climax  station  neutron  monitor  data  from 
the  University  of  Chicago.  These  data  are  regularly  downloaded  to  LSU  over  the 
NSI/DECnet  network,  and  plotted  for  use  in  the  selection  of  quiet  and  active  time 
periods  and,  eventually,  to  aid  in  determining  the  solar  modulation  level. 

Daily  reports  of  solar  activity  are  obtained  from  the  NOAA  Space  Environment 
Laboratory  again  using  the  NSI/DECnet  network.  These  reports  which  include  X- 
ray  fluence  plots,  geomagnetic  storm  and  solar  activity  summaries  and  forecasts,  and 
listings  of  recent  solar  events  are  printed  and  archived.  In  addition,  the  event 
reports  are  stored  online  for  later  "real-time"  access.  Geomagnetic  and  Solar  Indices 
such  as  Kp  and  Dst  are  obtained  from  the  National  Geophysical  Data  Center  on 
floppy  disks.  This  are  read  on  a  IBM-PC  and  transferred  to  the  online  VAX  database 
using  the  local  DECnet  network.  These  indices  are  used  to  correlate  the  active 
particle  fluxes  measured  with  ONR-604  to  geomagnetic  and  solar  disturbances. 

Proton  and  electron  particle  data  is  also  obtained  from  the  geosynchronous 
GOES  satellite  from  the  National  Geophysical  Data  Center.  This  data  is  received  in 
the  form  of  5  minute  averages  on  floppy  disk  and  for  the  full  time  resolution  (3 
seconds)  on  magnetic  tape.  In  addition,  the  5  minute  averages  are  received  every 
month  and  are  usually  only  a  month  or  so  off  real  time,  while  the  full  time 
resolution  dataset  can  be  delayed  up  to  8  months.  Finally,  the  OMNI  interplanetary 
conditions  (e.g.  solar  wind  velocity,  interplanetary  field)  are  obtained  from  the 
NSSDC  over  the  NSI/DECnet  network. 

It  is  expected  that  this  database  will  be  expanded  in  the  future  to  include  data 
from  other  CRRES  instruments  as  these  become  available  in  the  Science  Summary 
Database  maintained  online  at  AFGL.  In  particular,  proton  and  electron  data  may  be 
used  for  direct  comparison  with  the  IMP-8  and  GOES  data  and  for  completing  the 
low  charge,  low  energy  portion  of  the  GCR  Model. 
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5.  THE  HEAVY  ION  ENVIRONMENT  IN  LATE  1990-EARLY  1991 

a.  Overview 

The  CRRES  spacecraft  was  launched  during  a  period  near  the  maximum  of  the 
solar  activity  cycle,  as  illustrated  in  Figure  27  (top),  which  shows  the  last  four  solar 
cycles  as  seen  by  ground  based  neutron  monitors.  Data  from  both  Climax  and 
McMurdo  are  shown,  and  these  two  stations  "track"  well.  McMurdo,  located  near 
the  south  pole,  has  a  lower  geomagnetic  cutoff  than  Climax,  Colorado  but  the  mean 
energy  of  the  incident  particles  producing  the  counting  rate  is  almost  the  same. 

Note  that  90-91  was  a  period  in  which  the  neutron  monitors  reached  their  lowest 
level,  compared  to  previous  solar  cycles.  This  low  level  indicates  very  strong  solar 
modulation  for  the  cosmic  rays. 

The  bottom  portion  of  Figure  27  shows  an  expansion  of  the  late  1990-early  1991 
period.  CRRES  was  launched  just  after  the  minimum  in  the  neutron  monitor  rate, 
and  was  on  an  increasing  intensity  curve  (recovering  towards  solar  minimum 
conditions)  until  late  March  1991  when  the  flare  activity  halted  the  recovery. 
However,  the  time  period  between  ~1  Sept.  90  and  1  March  91  appears  to  be  a 
relatively  quiet  period  and  was  selected  as  the  first  six  month  period  for  model 
analysis. 

Figure  28  shows  the  selected  period  in  more  detail.  Compared  are  the  P2  and  P3 
rates  for  Apogee  (L  >  6.5)  and  for  all  proton  mode  data  (PM,  at  bottom)  and  the  GOES 
proton  counting  rates  P5  (39-82  MeV)  and  P3  (8.7-14.5  MeV).  Note  the  ONR-604  P2 
(>He)  rate  which  appears  relatively  constant.  No  major  heavy  ion  flares  occurred 
during  this  period.  The  ONR-604  P3  (L  >  6.5)  rate,  however,  does  show  some  proton 
or  electron  events.  These  are  correlated,  for  the  most  part,  with  the  low  energy 
proton  increases  seen  in  the  GOES  P3  channel.  The  larger  events  are  discernable  in 
the  GOES  P5  rate  as  well.  Now,  when  the  ONR-604  rates  are  plotted  for  all  proton 
mode  data  (bottom  two  plots),  the  P3  rate  shows  considerable  activity  while  P2 
remains  relatively  constant.  The  extra  spikes  in  P3  are,  presumably,  of 
magnetospheric  origin.  Only  near  apogee  can  the  ONR-604  counting  rate  for  P3  be 
used  to  trace  the  interplanetary  proton  component. 

The  top  panel  of  Figure  28  shows  the  interplanetary  quiet  times  (unshaded)  as 
determined  from  IMP-8  proton  data.  The  correlation  of  the  shaded  areas  with  the 
increases  seen  in  the  GOES  P3  rate  is  reasonably  good.  A  full  picture  of  the  IMP-8 
quiet  time  selections  for  1990  and  the  first  part  of  1991  is  shown  in  Figure  29,  where 
the  energies  corresponding  to  the  three  proton  and  three  helium  channels  are 
indicated  at  the  right.  Plotted  are  the  raw  counting  rates,  overlaid  with  the  quiet 
time  selections.  As  expected,  the  lower  the  energy  interval,  the  more  "active" 
appears  the  interplanetary  medium  and  the  smaller  are  the  intervals  selected  as 
"Quiet  Times."  For  heavy  ion  analysis,  however,  the  25-93  MeV/nudeon  Helium 
interval  is  most  applicable.  For  this  rate  there  are  only  a  few  "non-quiet"  periods, 
which  agrees  with  the  CRRES  F2  rate  ir  Figure  2P.  Thus,  for  the  selected  period 
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Figure  27.  (Top)  Four  solar  cycles  of  neutron  monitor  data  from  Climax  (courtesy 
University  of  Chicago)  and  McMurdo  (courtesy  Bartol  Research 
Institute).  (Bottom)  Expansion  of  the  latter  half  of  1990  and  the  first  part 
of  1991. 
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essentially  all  of  the  ONR-604  P2  and  P3  data  from  proton  mode  periods  can  be 
employed  for  environment  analysis. 
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Figure  29.  IMP-8  quiet  time  (clear  areas)  determined  from  six  different  channels. 

Figure  30  (top)  shows  the  number  of  PI  events  recorded  as  a  function  of  orbit 
and  the  L  distribution  of  PI  events  (bottom).  The  22-24  March  flare  period  starts  at 
orbit  584,  for  which  no  noticeable  increase  is  observed.  There  is  a  gradual  increase  in 
the  overall  rate  which  is  consistent  with  the  decreasing  modulation  shown  by  the 
neutron  monitors  (Figure  27).  In  L  Space,  most  of  the  events  occur  at  L  >  5  with  a 
band  at  L  =  6.5  which  is  the  location  of  the  geomagnetic  equator.  The  CRRES  orbital 
plane  and  the  geomagnetic  equatorial  plane  are  effectively  aligned  at  -orbit  500. 

What  is  interesting  is  that  there  are  a  handful  of  events  observed  at  low  L  values. 
These  are  most  likely  high  energy  GCR  events  but  could  also  be  pseudo-trapped 
particles  or  sub-threshold  events,  as  have  been  reported  by  others  (Oschlies  et  al., 
1989;  Durgaprasad  et  al.,  1990;  Adams  et  al.,  1991). 

During  initial  processing  of  the  ONR-604  data,  each  event  is  assigned  a 
preliminary  charge  and  a  preliminary  energy.  The  preliminary  charge  and  energy 
distributions  are  shown  in  Figure  31.  The  Mg  shows  the  largest  abundance  as 
expected  since  part  of  the  Ne  is  in  priority  2.  The  iron  peak  is  also  evident  as  are  a 
few  events  of  Co  +  Ni.  Note,  however,  that  there  are  events  above  Ni  and  events 


52 


Counts 


140 


CRRES  Orbit  Number 


100  150  200  250  300  350  400  450  500  550  600  650  700 

CRRES  ORBIT  NUMBER 

Figure  30.  The  number  of  PI  events  as  a  function  of  orbit  number  (top)  and  the 
L  distribution  of  PI  events  (bottom). 
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below  Oxygen.  Also,  the  odd-Z  elements  are  more  abundant  than  expected.  This 
suggests  that  (1)  the  preliminary  ZCAL  charge  assignment  is  not  completely 
accurate,  and  (2)  there  is  potential  "background",  at  the  <15%  level,  among  the  PI 
events.  The  energy  distribution  shows  the  expected  energy  ranges  (c.f.  Figure  9)  for 
each  element,  with  the  exception  of  the  low-Z  events  which  should  not  appear  in 
PI.  There  are  only  a  few  events  outside  the  nominal  energy  limits. 

The  preliminary  Z  values  are  calculated  from  only  the  information  in  the  last 
two  detectors  triggered.  Thus,  if  a  particle  interacts  in  the  stack  sending  a  low-Z 
fragment  forward,  it  is  the  fragment’s  Z  that  will  be  assigned.  This  can  account  for 
some  of  the  low-Z  events,  but  probably  not  all.  Multiple  particles  are  also  a 
possibility  as  are  bit  errors  that  assign  a  P2  event  to  the  PI  category.  Considerably 
more  investigation  is  required  both  to  develop  a  better  charge  calibration  and  to 
understand  the  background  contributions  among  the  PI  events. 


b.  Solar  Modulation 

It  is  well  known  (see  reviews  by  Jokipii  1971,  Fisk  1979,  McKibben,  1986)  that  at 
energies  below  a  few  GeV  per  nucleon  the  intensity  and  the  spectral  shape  of  the 
cosmic  rays  arriving  at  the  orbit  of  Earth  are  significantly  modified  by  solar 
modulation.  The  cosmic  rays  diffuse  into  the  heliosphere  against  the  outward¬ 
flowing  solar  wind,  carrying  frozen-in  interplanetary  magnetic  field,  which  convects 
the  cosmic  rays  out  of  the  heliosphere.  In  this  process,  the  cosmic  rays  lose  energy  to 
the  expanding  field  (adiabatic  deceleration),  and  their  energy  spectrum  is  modified. 

A  spherically  symmetric  model  of  solar  modulation  has  been  developed 
(Parker  1965;  Jokipii  1971;  Urch  and  Gleeson  1972;  Fisk  1979)  which  explains  most  of 
the  gross  features  of  the  modulation  process.  This  model  includes  the  effects  of 
diffusion,  convection,  and  adiabatic  deceleration  (but  not  drifts  due  to  the  gradient 
and  curvature  of  the  magnetic  field)  and  assumes  that  these  three  physical  processes 
are  in  equilibrium  in  the  heliosphere.  Quantitatively  these  effects  are  represented 
by  a  Fokker-Planck  equation  in  which  the  parameters  are  the  solar  wind  velocity, 
the  diffusion  coefficient  and  the  radius  of  the  heliosphere,  with  the  cosmic-ray 
differential  energy  spectrum  in  local  interstellar  space  as  a  boundary  condition. 

Evenson  et  al.  (1983)  (see  also,  Garcia-Munoz  et  al.  1986),  solving  this  equation 
numerically,  have  analyzed  the  simultaneous  modulation  of  electrons,  protons,  and 
helium  nuclei  over  the  1965-1979  period  involving  more  than  one  solar  cycle.  They 
find  that  in  general  the  model  fits  the  data  quite  well.  In  this  model  the  degree  of 
modulation  at  a  heliospheric  radius  r  is  given  by  the  modulation  parameter 
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where  V(r')  is  the  solar  wind  velocity,  K(r')  is  the  radial  part  of  the  diffusion 
coefficient,  and  R  is  the  radius  of  the  heliosphere.  An  insight  into  the  physical 
meaning  of  4>  is  obtained  from  the  "force-field"  approximation  (Gleeson  and  Axforo 
1968)  in  which  O  corresponds  to  a  "potential  energy,"  that  in  the  particular  case  in 
which  the  diffusion  coefficient  is  proportional  to  particle  rigidity  takes  the  simple 
form 


<D=|ZeU(r),  (2) 

where  Ze  is  the  particle  charge.  This  potential  energy  has  been  identified  as  the 
mean  energy  loss  that  the  particles  experience  in  penetrating  the  heliosphere  to  a 
radius  r.  The  parameter  may  be  expressed  in  rigidity  or  energy  units,  connected 
through. 


O  (MV)  =  (MeV/ nucleon) 

The  modulated  nucleonic  differential  energy  spectra  obtained  by  the  numerical 
solution  of  the  Fokker-Planck  equation  are  almost  completely  determined  by  the 
value  of  the  modulation  parameter  <j>.  Combinations  of  the  parameters  R,  r,  V(r), 
and  K(r)  giving  the  same  value  of  $  will  lead  to  modulated  nucleonic  spectra  which 
are  very  nearly  equal  to  each  other  (Urch  and  Gleeson  1972).  Therefore,  values  of 
the  modulation  parameter  $  are  used  to  specify  different  levels  of  modulation,  as 
was  shown  in  Figure  4  for  calculated  Hydrogen  and  Helium  spectra  for  different 
values  of  the  parameter  O  . 

Figure  32  (from  Garcia-Munoz  et  ah,  1987)  shows  the  time  dependence  of  the 
modulation  parameter  used  by  Evenson  et  al.  (1983)  for  1973-1980,  compared  with 
the  time  dependence  of  the  Climax  Neutron  Monitor  count  rate  (given  as  the 
percent  decrease  below  the  1954  solar  minimum  level).  The  period  1974-1978, 
corresponding  to  solar  minimum  conditions,  was  characterized  by  a  nearly  constant 
level  of  solar  modulation  for  which  the  average  value  of  the  modulation  parameter 
was  about  <|>  =  490  MV,  or  <1>  =  245  MeV  per  nucleon  for  A/Z  =  2  particles.  This  time 
period,  during  which  much  cosmic  ray  data  was  recorded,  becomes  a  baseline  to 
which  more  recent  data,  in  particular  the  CRRES  results,  are  compared. 

There  are  several  approaches  to  determining  the  solar  modulation  level  for  a 
specific  time  period  during  the  CRRES  mission: 

(a)  Demodulation  of  CRRES  elemental  spectra 

(b)  Demodulation  of  IMP-8  elemental  spectra 

(c)  CRRES  Heavy  Ion  Rate  Modeling 

(d)  Neutron  Monitor  Correlation 
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Figure  32.  The  modulation  parameter  <t>  (MV)  for  1974-1980  determined  by 
Evenson  et  al.  (1983)  from  the  flux  of  70-95  MeV/ nucleon  Helium 
compared  to  (B)  the  Climax  Neutron  Monitor  rate. 

The  normalization  procedures  and  background  subtraction  methods  for  ONR-604 
have  not  yet  been  completed,  so  that  actual  flux  values,  which  could  be  fit  by  the 
solar  modulation  model  described  above,  are  not  available.  However,  IMP-8  can  be 
used  in  the  same  way,  and  some  preliminary  IMP-8  data  has  been  provided  by 
Chicago.  Technique  (d)  exploits  the  correlation  shown  in  Figure  32,  developed  for 
the  previous  solar  cycle.  It  is  not  known  whether  the  same  correlation  holds  for  the 
present  cycle,  and  this  approach  remains  to  be  investigated.  What  we  have 
concentrated  upon  during  the  past  year  is  method  (c),  which  will  be  described  below, 
and  using  (b)  for  validating  or  comparing  the  results  obtained. 

The  ONR-604  PI  and  P2  rates  are  dominated  by  heavy  ions  —  Mg  -  Fe  for  PI  and 
He  -  O  for  P2.  Both  rates  respond  mainly  to  GCR  and  are  relatively  insensitive  to 
magnetospheric  populations,  at  least  during  the  apogee  portion  of  the  orbit  as  was 
indicated  by  Figure  28.  Thus,  using  the  instrument  description  and  operating  modes 
described  in  section  2,  it  should  be  possible  to  calculate  the  expected  priority  rates  for 
any  given  particle  spectrum  outside  the  spacecraft. 

What  has  been  done  is  the  following.  For  each  element,  the  energy  range 
appropriate  to  a  given  logic  condition  is  determined  (c.f.  Figure  9).  Then  the 
geometry  factor  is  calculated  as  a  function  of  energy  (range)  as  illustrated  in  the  left 
part  of  Figure  33.  Any  incident  spectrum  can  then  be  integrated  over  the 
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Figure  33.  Energy  dependent  geometry  factors  for  Ne  and  Fe  events  (left)  and  the  geomagnetic  transmission 
function  (right). 


appropriate  energy  range,  folding  in  the  geometry  factor,  to  predict  the  number  of 
events  per  second.  Repeating  the  process  for  all  of  the  elements  in  a  particular  logic 
condition  (such  as  PI),  and  summing  the  results,  gives  a  calculated  rata 
(counts /second)  to  compare  to  the  flight  data. 

The  input  spectra  for  such  a  calculation  are  the  differential  energy  spectra 
outside  the  instrument,  i.e.  at  the  CRRES  orbit.  This  is  obtained  by  folding  the 
interplanetary  spectra  with  the  geomagnetic  transmission  function  (see  the  right 
portion  of  Figure  33).  For  data  limited  to  the  apogee  portion  of  the  orbit  (proton 
mode  for  ONR-604)  there  is  little  geomagnetic  inhibition  except  at  the  lowest 
particle  energies.  The  interplanetary  spectra  are  calculated  by  modulating  the  spectra 

in  local  interstellar  space  (LIS)  to  the  orbit  of  Earth.  In  this  process  d>  is  the 
parameter,  so  that  the  calculations  can  be  performed  repeatedly  for  a  variety  of 
values  of  <I>.  The  value  of  O  that  gives  calculated  rates  in  accord  with  the  data  is 
then  the  appropriate  level  of  modulation  for  the  time  interval  being  studied. 

Figure  34  shows  an  example  of  such  a  study  for  PI  and  P2  rates.  Plotted  are 

calculated  counting  rates  for:  PI  =  PIS,  PI*  =  PISA,  P2  =  P2AS,  P2*  =  P2AS  (Z  >  2), 
which  correspond  to  particles  not  entering  or  leaving  the  sides  of  the  experiment 
and,  except  for  PI,  stopping  in  the  detector  stack.  The  curves  are  the  result  of  many 
calculations,  each  for  a  different  level  of  modulation,  shown  on  the  horizontal  axis 
in  units  of  MeV/nucleon.  Solar  minimum  conditions  exhibit  a  residual 
modulation  level  of  ~250  MeV/nucleon,  so  the  lower  values  of  O  will  not  be 
applicable,  while  previous  solar  maximum  conditions  have  given  O  values  of  600- 
800  MeV/nucleon.  The  PI*  and  P2*  rates  respond  principally  to  low  energy 
(stopping)  Mg-Fe  and  C-O  nuclei,  respectively,  and  are  slightly  more  sensitive  to 

changes  in  <I>  than  is  PI  which  includes  penetrating  particles  as  well.  The  P2  rate  is 
shown  for  comparison  but  is  not  used  since  the  helium  region  of  P2  events  shows 
some  background  contamination. 

The  equivalent  counting  rates  must  be  obtained  from  the  pulse  height  analyzed 
events,  since  the  hardware  logic  circuitry  does  not  conform  to  any  of  the  conditions 
specified  above.  Each  event  is  read  from  tape  and  its  discriminator  bit  pattern  is 
compared  to  the  conditions  described  previously.  In  the  case  of  P2*,  a  preliminary 
charge  assignment  of  ZCAL  >  2.8  is  required  as  well.  Each  event  meeting  the  criteria 
is  counted  over  a  specified  active  time  interval  to  obtain  the  "measured"  counting 
rate.  These  are  shown,  for  PI,  PI*  and  P2*,  at  the  right  side  of  Figure  34.  Note  that 
the  three  rates,  when  extended  to  the  appropriate  curve,  all  give  a  relatively 
consistent  value  for  the  solar  modulation  parameter,  0=  550  ±  50  MeV/nucleon. 

While  Figure  34  shows  a  consistent  result  for  <J>,  it  remains  to  be  demonstrated 
that  this  value  is  correct.  Such  a  validation  can  be  provided  by  IMP-8  data  (in  the 
absence  of  normalized  CRRES  spectra)  since  the  IMP-8  instrument  is  well 
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Calculated  priority  counting  rates  for  ONR-604  (left)  and  comparison  of  the  calculations 


understood  and  has  been  utilized  for  many  years  in  such  studies.  Figure  35  (top) 
shows  IMP-8  oxygen  spectra  provided  by  M.  G.  Munoz  for  the  late  199C-early  1991 
period,  as  well  as  for  previous  periods,  to  be  used  in  comparative  studies.  Beginning 
with  the  PLEXP  model  for  local  interstellar  spectra,  the  lower  set  of  curves  in  Figure 
35  (bottom)  shows  the  solar  modulation  prediction  compared  to  the  IMP-8  data.  A 
best  fit  value  of  650  ±  30  MeV/nucleon  (1300  ±  60  MV)  is  obtained,  somewhat  above 
the  values  for  C>  inferred  from  the  priority  rate  analysis.  Assuming  that  the  IMP-8 
level  is  correct,  then  there  must  be  a  problem  with  the  rate  model  or  with  the  ONR- 
604  data  itself.  Alternatively,  the  PLEXP  local  interstellar  spectra  may  not  be  the 
most  applicable,  and  this  would  affect  the  curves  shown  in  both  Figures  34  and  35 
(bottom). 

There  are  other  calculations  of  the  local  interstellar  spectra,  each  corresponding 
to  a  different  pathlength  distribution  for  galactic  propagation  (see  the  discussion  in 
the  next  section).  Of  course,  each  model  must  reproduce  the  solar  minimum  spectra 
that  were  shown  in  Figure  3,  so  the  differences  in  LIS  are  small.  The  effect  of  these 
different  local  interstellar  spectra  on  the  calculated  PI,  PI*  and  P2*  counting  rates  is 
shown  in  Figure  36  where  each  of  the  curves  corresponds  to  a  different  model.  Note 
that  the  effect  on  the  calculated  rates  in  the  <t>  region  of  interest  is  very  small, 
corresponding  to  introducing  an  uncertainty  of  ±50  MeV/nucleon  into  the  inferred 
value  of  0. 

The  geomagnetic  transmission  function  has  been  checked,  as  has  the  energy 
dependent  acceptance,  and  no  major  uncertainties  have  been  found.  Since  the  rate 

analysis  shows  a  lower  value  of  <J>,  the  measured  counting  rate  should  be  lower  to 
bring  the  inferred  <b  values  into  agreement.  This  might  indicate  a  residual 
background  among  the  PI  and  P2  events  which  has  not  been  removed  by  the  strict 
logic  cuts  employed.  Further  investigation  will  be  required  to  try  to  resolve  the 
apparent  discrepancy.  However,  for  the  purposes  of  this  report  a  preliminary  value 
of  =  625  ±  75  MeV/nucleon  is  adopted. 


c.  Galactic  Cosmic  Rays 

The  Galactic  Cosmic  Rays  (GCR)  are  a  continuous  source  of  ionizing  radiation 
of  extra-solar  origin.  For  the  past  several  decades  the  composition  and  energy 
spectrum  of  these  high  energy  nuclei  have  been  measured  by  a  number  of  satellite 
experiments  and  ONR-604  represents  the  latest,  current  technology  instrument  in 
this  series.  Thus,  the  GCR  Model  to  be  developed  under  this  effort  has  a  wealth  of 
previous  measurements  which  can  be  used  as  a  baseline.  In  fact,  ONR-604  will  not 
provide  measurements  over  the  full  charge  and  energy  range  required  by  the  Model, 
so  the  baseline  dataset  will  be  needed  to  set  and/or  constrain  a  majority  of  the  Model 
parameters. 
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Figure  35.  Oxygen  energy  spectra  from  IMP-8  for  different  time  periods  (top)  and 
*  modulation  fits  to  the  IMP=8  data  for  1978  and  1990-91  (bottom). 


Figure  36.  Comparison  of  calculated  PI,  PI*  and  P2*  rates  for  different  local  interstellar  spectra. 


The  current  GCR  baseline  database  contains  close  to  1000  data  points  from  mor; 
than  50  literature  references  and  covers  the  most  important  elemental  ratios  and 
spectra.  The  database  was  initially  developed  in  the  early  1980's  and  a  majority  of  th 
measurements  available  then  were  made  during  the  extended  solar  minimum 
period  of  1972  - 1978.  Since  then  relatively  few,  new  missions  have  been  flown  and 
the  statistics  collected  during  the  current  solar  cycle  have  been  hampered  by  the  two 
relatively  deep  solar  maximums  bracketing  a  rather  narrow  solar  minimum  period 
(c.f.  Figure  27). 

Recently,  some  new  data  has  become  available  (e.g.  Seo  et  al.,  1991;  Engelmann 
et  al.,  1990;  Ferrando  et  al.,  1991),  and  these  need  to  be  incorporated  into  the  database. 
However,  these  measurements  were  made  under  solar  modulation  conditions,  very 
different,  in  some  cases,  from  the  measurements  in  the  baseline  database.  For 
example,  Ferrando  et  al.  (1991)  r .  rts  on  low  energy  (-100  MeV/nudeon)  GCR 
composition  measurements  made,  uring  1986-1987  with  the  High  Energy  Telescope 
onboard  the  Voyager  2  spacecraft  /hile  the  spacecraft  was  at  22  AU.  While  the  time 
period  involved  corresponds  to  the  last  solar  minimum,  these  outer  heliospheric 
measurements  can  not  be  directly  compared  with  those  made  at  the  orbit  of  Earth. 
Thus,  we  will  need  to  devote  some  effort  to  evaluating  new  measurements  before 
including  them  in  the  database. 

A  portion  of  the  energy  spectra  measurements  from  Figure  3  are  shown  in 
Figure  37,  for  He,  C  and  Fe  data  with  curves  from  initial  GCR  Model  calculations 
These  results  used  a  calculated  local  interstellar  spectrum  (LIS)  modulated  by  the 
technique  described  earlier  with  solar  minimum  conditions  corresponding  to  the 
1972  - 1978  period  (<b  =  490  MV).  The  LIS  was  calculated  using  a  computer  program 
that  is  based  upon  the  weighted-slab  GCR  interstellar  propagation  technique  (Fichtel 
and  Reames,  1968)  and  which  is  fully  described  in  Garda-Munoz  et  al.  (1987).  In 
essence,  this  technique  starts  with  a  set  of  relative  abundances  of  elements  and 
isotopes  at  the  GCR  source,  assumes  a  common  source  energy  spectrum,  as  well  as  a 
particular  interstellar  medium  (ISM)  composition  and  density,  and  then  calculates 
new  abundances  and  spectra  for  a  series  of  discrete  ISM  slabs.  This  calculation 
includes  the  effects  of  production  and  loss  of  spedes  via  nuclear  interaction, 

radioactive  decay  (P+,  P‘,  and  e*  capture),  electron  capture  and  loss,  and  ionization 
energy  loss.  These  "slab  results"  are  then  integrated  over  an  assumed  pathlength 
distribution  (PLD)  to  obtain  the  LIS. 

One  might  note  at  this  point  that  the  "Slab  Calculations"  involve  only  atomic 
and  nuclear  effects  which  can,  in  prindpal,  be  measured  in  the  laboratory,  while  the 
PLD  involves  the  astrophysical  details  of  GCR  interstellar  propagation  and  must  be 
either  theoretically  derived  or  obtained  from  the  GCR  data.  Thus,  obtaining  a  good 
representation  of  the  Galadic  Cosmic  Rays  is  actually  a  multi-parameter  problem 
and  it  is  important  to  use  laboratory  measurements,  GCR  data  and  other 
astrophysical  information  to  constrain  as  many  of  these  parameters  as  possible.  For 
example,  of  particular  importance  are  the  nuclear  interaction  cross  sections,  but  only 
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Figure  37.  Measurements  of  GCR  He,  C  and  Fe  flux  from  the  baseline  database 
along  with  predictions  (solid  curve)  of  the  initial  GCR  Model. 

a  few  interaction  channels  have  been  measured  over  the  energy  range  required  by 
the  cosmic  ray  problem  (Guzik,  1990).  Therefore,  there  is  a  heavy  reliance  on 
calculated  cross  sections  (e.g.  Silberberg  and  Tsao,  1990)  but  these  can  differ  from 
measured  values  by  up  to  60%  (Webber,  Kish  and  Schrier,  1990).  Constraining  the 
Model  parameters  is,  thus,  a  complex  process  involving  direct  laboratory 
measurements,  semi-emperical  calculations  and  good  judgement.  Considerable 
effort  was  required  to  establish  the  parameter  database  for  the  initial  GCR  Model 
calculations  and  additional  effort  will  be  required  to  maintain  this  database  as  new 
laboratory  measurements  or  other  constraints  become  available. 

The  PLD  used  for  the  results  shown  in  Figure  37  is  an  exponential  distribution 
of  pathlengths  with  a  mean  that  is  a  decreasing  power  law  in  total  energy.  This 
distribution,  referred  to  as  the  Power-Law  Exponential  (PLEXP),  is  consistent  with  a 
homogeneous  model  of  cosmic  ray  propagation  which  allows  a  small,  but  increasing 
"leakage"  from  the  "galactic"  confinement  volume  with  increasing  energy.  These 
results  appear  to  fit  the  data  of  Figure  37  quite  well.  Helium,  Carbon  and  Iron, 
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however,  are  "primary"  species.  That  is,  their  relative  abundance  in  the  measured 
GCR  is  dominated  by  their  abundance  at  the  source  and  is  not  significantly  affected 
by  the  "secondary"  component  generated  by  the  nuclear  fragmentation  of  heavier 
species.  Therefore,  primary  species  reflect  conditions  at  the  cosmic  ray  source  rather 
than  conditions  during  GCR  propagation.  Pure  secondary  species,  which  have  a 
negligible  source  component,  are  very  sensitive  to  the  exact  form  of  the  PLD.  In  fact, 
by  examining  ratios  of  secondary  to  primary  elements  the  PLD  parameter  in  the 
GCR  Model  can  be  constrained. 

Figure  38  shows  measurements  from  the  GCR  database  for  the  B/C  secondary 
to  primary  ratio,  along  with  the  GCR  Model  results  using  the  PLEXP  PLD  (long 
dashed  curve).  For  energies  greater  than  1  GeV/nudeon  the  prediction  represents 
the  data,  but  overestimates  the  measurements  for  lower  energies.  Thus,  while  the 
Model  may  correctly  predict  the  primary  spectra  it  will  overpredict  the  flux  of 
secondary  GCR  at  100  MeV/nucleon  by  -60%.  An  alternative  PLD  referred  to  as  the 
Double-Sided  Power  Law  (DSPL)  pathlength  distribution  is  shown  as  the  solid  curve 
in  Figure  38  and  provides  a  much  better  fit  to  the  ratio.  The  DSPL  has,  in  essence, 
the  same  form  as  the  PLEXP  for  high  energy  (>  1  GeV/nudeon),  but  for  low  energy 
the  exponential  mean  is  a  decreasing  function  of  decreasing  energy.  Such  a  PLD 
may  be  consistent  with  propagation  models  which  involve  a  "galactic  wind"  that 
flows  outward  from  the  galactic  plane  in  a  fashion  similar  to  the  solar  wind  which 
flows  outward  from  the  Sun  (Garcia-Munoz  et  al.,  1987). 

Further  details  of  the  PLD  can  be  constrained  by  requiring  multiple  secondary 
to  primary  ratios  widely  separated  in  charge  to  be  simultaneously  fit.  Figure  39 
shows  the  collection  of  data  for  the  ratio  Sub-Fe/Fe  where  the  "Sub-Fe"  group 
includes  the  elements  Sc,  Ti,  V,  Cr  and  Mn.  The  DSPL  results  which  fit  the  B/C 
ratio  is  shown  as  the  solid  curve  and  is  quite  likely  low  of  the  data  for  energies  less 
than  -1.0  GeV/nudeon.  This  discrepancy  has  been  used  to  argue  that  the  PLD  must 
be  depleted  of  short  pathlengths  (Garda-Munoz  et  al.,  1987).  Since  the  Iron  group 
species  have  a  larger  total  nuclear  interaction  cross  section,  and  hence  a  shorter 
mean  interaction  length,  than  the  Carbon  group  nudei,  the  iron  secondaries  will  be 
more  sensitive  to  the  short  pathlength  characteristics  of  the  PLD.  By  removing  or 
"truncating"  the  short  pathlengths  from  the  PLD  the  calculation  will  cause  more 
Iron  nuclei  to  interact  and  raise  the  Sub-Fe/Fe  ratio  (dashed  curves  in  Figure  39) 
while  not  affecting  the  prediction  for  B/C  (short  dashed  curve  in  Figure  38). 

The  spectra  calculated  using  these  various  PLDs  are  shown  in  Figure  40.  The 
solid  curves  are  the  same  as  in  Figure  37,  that  is  using  the  PLEXP,  the  short  dashed 
curves  are  results  using  the  DSPL  and  the  long  dashed  curves  result  from  a 
"truncated"  PLD.  For  Helium  and  Carbon  the  three  PLDs  yield  almost  identical 
results,  but  the  differences  are  apparent  in  the  low  energy  (<700  MeV/nucleon)  Iron 
spectrum.  As  energy  deposit  is  proportional  to  Z2/E  an  error  made  in  calculating 
the  low  energy  Iron  spectrum  may  translate  into  a  relatively  large  uncertainty  in  the 
Linear  Energy  Transfer  (LET)  spectrum. 
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38.  Measurements  of  B/C  with  Model  predictions  using  different  PLDs. 
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Figure  39.  Measurements  of  Sub-Fe/Fe  with  Model  predictions  using  different  PLDs 
(see  text). 
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Figure  40.  Spectra  predicted  by  the  GCR  Model  for  different  PLDs  (see  text). 

In  fact,  the  accuracy  to  which  the  PLD  can  be  determined  and,  consequently,  the 
spectra  can  be  calculated,  depends  upon  the  accuracy  of  the  other  Model  parameters. 
For  example,  while  the  solid  curve  in  Figure  39  (no  truncation)  appears  to 
underpredict  the  low  energy  data,  it  can  be  argued  that  the  dashed  curves  (with 
truncation)  are  high  for  all  energies.  This  could  be  due  to  having  an  incorrect 
source  abundance  for  the  Sub-Fe  species  in  the  Model,  and  reducing  these 
abundances  would  lower  the  dashed  curves.  In  addition,  PLD  truncation  is  still 
controversial  and  the  case  for  truncation  is  based  largely  on  the  low  energy  IMP-8 
measurement  (filled  diamond  in  Figure  39).  The  recent  measurement  by  Ferrando 
et  al.  (1991)  provides  the  only  other  data  in  the  IMP-8  charge  and  energy  range. 

Their  value  for  Sub-Fe/Fe  is  about  30%  less  than  the  IMP-8  measurement  and  is  in 
good  agreement  with  the  solid  curve  (no  truncation)  of  Figure  39.  As  discussed 
above,  however,  it  is  not  immediately  clear  whether  the  Ferrando  et  al. 
measurement  made  with  an  instrument  onboard  the  Voyager  2  spacecraft  at  22  AU 
can  be  so  directly  compared  with  the  IMP-8  data  made  in  the  vicinity  of  Earth. 
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Finally,  the  nudear  interaction  cross  sections  for  the  production  of  Sub-Fe  species 
from  Iron  are  also  largely  unmeasured  and,  thus,  the  Model  uses  mostly  calculated 
cross  sections.  If  these  calculated  cross  sections  are  increased  in  value  then  the  solid 
curve  in  Figure  39  would  move  up  and  it  may  be  possible  to  fit  the  Sub-Fe/Fe  and 
B/C  data  simultaneously  without  requiring  a  truncated  pathlength  distribution. 

Refining  the  GCR  Model  and  evaluating  the  consequences  of  these 
uncertainties  is  a  time  consuming  process  involving  multiple  calculations.  First  an 
alternate  set  of  input  parameters,  such  as  the  source  abundances,  source  spectra  or 
nuclear  cross  sections,  must  be  determined  from  the  literature,  derived  or 
calculated.  Next  the  "Slab  Calculation"  program,  which  takes  roughly  10  to  12  CPU 
hours  on  an  IBM  3090,  must  be  run  with  the  new  inputs.  Using  these  new  slab 
results  the  PLD  parameters  must  be  adjusted  so  that  the  B/C  and  Sub-Fe/Fe  data  are 
simultaneously  fit.  The  final  results  of  this  process  are  then  compared  with  the 
spectra  and  ratio  measurements  from  the  baseline  database.  If  further  adjustments 
are  still  called  for  then  the  process  begins  again.  Usually  3  to  4  iterations  through 
this  loop  and  several  months  of  man-effort  are  needed  to  achieve  an  additional 
level  of  refinement  in  the  GCR  Model. 


d.  Preliminary  Predictions 

The  preceding  sections  have  discussed  the  needed  inputs  to  the  Interplanetary 
Heavy  Ion  model  and  the  current  state  of  the  analysis.  While  it  is  still  premature  to 
issue  a  definitive  environment  model,  a  preliminary  "working"  prediction  can  be 
made.  Since  the  time  period  of  interest  (1  Sept.  90  -  28  Feb.  91)  contains  no  large 
solar  flares,  a  flare  component  need  not  be  considered.  Thus,  the  first  period 
environment  is  dominated  by  the  Galactic  Cosmic  Rays. 

Using  the  PLEXP  spectra  in  Local  Interstellar  Space  and  a  modulation  level  of  <I> 
=  625  MeV/ nucleon,  the  LPARL  group  (D.  L.  Chenette,  lead  investigatory)  has 
calculated  the  flux  of  cosmic  rays  at  the  orbit  of  Earth  given  in  Table  2  and  illustrated 
in  Figure  41  for  the  elements  H,  He,  O  and  Fe. 

To  investigate  the  effect  of  these  particles  on  components,  it  is  necessary  to 
calculate  the  Linear  Energy  Transfer  (LET)  distribution  for  different  LET  thresholds 
and  behind  different  levels  of  shielding.  This  has  been  done  for  shielding 
thicknesses  of  50,  100,  200,  400,  800  and  1600  mils  of  Aluminum  for  a  grid  of  different 
thresholds.  The  numerical  results  are  given  in  Table  3  and  are  shown  graphically  in 
Figure  42.  These  preliminary  results  can  be  utilized  by  the  Product  Associated 
Working  Groups. 

What  has  not  been  completed  is  a  study  of  the  variation  in  spectra  or  LET 
distributions  introduced  by  the  uncertainty  in  the  level  of  solar  modulation  or  the 
different  US  spectra  that  were  described  above.  Such  a  study  will  indicate  the  range 
of  possible  variation  about  the  "working"  values  in  Tables  2  and  3,  and  will  allow 


69 


Differential  flux  (m2  s  sr  MeV)-i 


1.E+01  1.E+02  1.E+03  1.E+04  1.E+05 

Kinetic  Energy  (MeV/nucleon) 

Figure  41.  Predicted  energy  spectra  for  GCR  nuclei  H,  He,  O  and  Fe  for  the  late  ’90- 
early  VI  time  period. 
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Figure  42.  Integral  number  of  particles  with  LET  greater  than  threshold,  LO,  as  a 
-  function  of  LO.  The  curves  are  labeled  with  the  thickness  of  Aluminum 

shielding  (in  mils)  assumed  in  the  calculation. 
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the  "user"  of  these  results  to  estimate  the  uncertainties  in  his  analysis.  The  study  of 
these  variations  will  be  one  of  the  next  tasks  undertaken  by  the  modeling  group. 


6.  SUMMARY 

The  CRRES  mission  provides  an  opportunity  for  new,  detailed  in  situ 
observations  of  the  Geospace  radiation  environment  and  its  effects  on 
microelectronic  components.  The  ONR-604  experiment  provides  the  monitoring  of 
high  energy  heavy  ions  and,  together  with  previous  and  contemporary  data, 
provides  a  means  for  describing  the  radiation  environment  experienced  by  the 
CRRES  spacecraft.  This  environment  is  being  codified  and  described  by  an 
Interplanetary  Heavy  Ion  model  that  has  been  subject  of  this  interim  technical 
report. 

The  overall  modeling  effort  as  well  as  the  operation  and  performance  of  the 
ONR-604  instrument  through  pre-launch  calibrations  and  initial  on-orbit 
operations  has  been  described  in  detail  as  necessary  background  information  for  the 
data  analysis  and  modeling  effort.  The  data  reduction  and  processing  system 
development  has  encompassed  the  reading,  checking,  formatting  and  splitting  of 
the  data  into  databases  as  well  as  the  capability  to  plot  and  correlate  the  data  in  a 
number  of  ways.  Auxiliary  databases  have  been  established  and  are  used  in  the 
modeling  effort.  What  remains  is  to  develop  the  background  subtraction  and 
normalization  portion  of  the  system,  in  order  to  convert  the  ONR-604  results  into 
meaningful  flux  values  as  a  function  of  charge  and  energy.  In  addition,  solar  flare 
analysis  and  modeling  routines  have  not  yet  been  finished. 

The  analysis  of  the  environment  data,  both  CRRES  and  other  data,  has 
proceeded  to  the  point  of  establishing  a  preliminary,  "working"  environment 
model  for  the  period  ~1  Sept.  90  - 1  March  91,  characterized  by  relatively  quiet 
interplanetary  conditions.  This  environment  is  characterized  by  a  high  level  of 

solar  modulation  (<J>  =  625  ±  75  MeV/nucleon)  and  correspondingly  depressed  levels 
of  Galactic  Cosmic  Rays. 

The  Solar  Modulation  process  and  the  techniques  developed  to  determine  the 
level  of  modulation  are  reviewed  in  detail,  revealing  that  the  two  techniques 
studied  give  similar  results,  qualitatively,  but  disagree  quantitatively  indicating  that 
further,  refined  analysis  is  required. 

Galactic  Cosmic  Ray  spectral  data  from  previous  epochs  are  summarized  and 
used  to  normalize  the  model  calculations.  New  data  is  described  as  well  as  the 
techniques  to  obtain  a  fully  self-consistent  set  of  spectra  in  Local  Interstellar  Space. 
Current  models  do  a  reasonably  good  job  in  reproducing  the  measured  data  but  are 
not  optimum  descriptions  of  the  environment.  Fur^er  refinements  are  needed  in 
the  GCR  model  to  bring  it  to  its  optimum  state  for  use  in  the  Interplanetary  Heavy 
Ion  Model  analysis. 
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The  "working"  environment  presented  here  in  terms  of  particle  Tux  and  LET 
spectra  can  be  used  by  the  community  as  a  preliminary  result  for  analyzing 
microelectronic  data.  It  is  hoped  that  the  product  associated  working  groups  will 
provide  both  analyses  and  suggestions  that  can  be  folded  into  the  process  described 
in  this  report  to  help  produce  a  better  product  for  the  user  community. 


* 
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Table  2:  Calculated  Galactic  Cosmic  Ray  Spectra  for  Z  =  1  to  28  and  E  =  50  to  over  60,000  MeV/nucleon 
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Table  2:  Calculated  Galactic  Cosmic  Ray  Spectra  for  Z=1  to  28  and  E  =  50  to  over  60,000  MeV/nucleon 
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Table  3:  Integral  LET  spectra  for  aluminum  shields  of  thickness  ranging  from 
50  mils  to  1600  mils  of  aluminum. 


6-Sep-9 1  Integral  LET  Spectra  for  August  1990  -  March  1991  CRRES  Cosmic  Ray  Model 
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